SHIP SYNTHESIS MODEL FOR 
COAST GUARD CUTTERS 



Michael Jay Goodwin 



» KifOX Ult^ABX 
post TB AOUATE SCHOOt 
•»0*IIWEY. CALIFORNIA 93940 



SKI? SYNTHESIS r.iODSL FOR 



COAST GUARD CUTTERS 

MIvKASL JAY GOOD^VIN 

/•" 

B.S., Webb Institute 

1969 

SUBMITTED IN PARTIAL FULFILLBEK 
OF THE RE QUIREI.EKTS FOR THE 
DEGREE OF OCEAN ENGINEER 
and 

TKS DEGREE OF MASTER OF 
SCIENCE 
at the 

I^iASSACHUSSTTS INSTITUTE CF 
TECHNOLOGY 
June, 1975 



SS^i* 




«ONT£,-!cY. C«i.lr 



-i.OLi, 

-J340 



ABSTRACT 



Ship Synthesis Model for Coast Guard Cutters 
byj Michael Jay Goodv/in 

Submitted to the Department of Ocean Engineering on 

May 9» 1975# in partial fulfillment of the requirements 

for the degree of Ocean Engineer and the degree of 

Master of Science in Naval Architecture and Marine Engineering 



In order to determine the optimum design of a Coast 
Guard cutter, it is necessary to perform a large number 
of feasibility studies. This is an impossible task if 
done using hand calculations. A mechanized version of 
these hand calculations is developed which will permit 
the calculations to bo p>erformod by a digital computer. 

The final progratm listing is given together with the 
development of the empirical relations used. 

The program is applicable to search and rescue and 
patrol type cutters with lengths between 150 and 400 feet. 
Twin screw propulsion is assumed. 

A check on the accuracy of the program is made using 
three U.S. Coast Guard cutters now in service and the 
U.S. Navy's Patrol Frigate design. Tho program is found to 
be much more accurate than existing programs v.’hen applied 
to Coast Guard vessels. The program predicts the full load 
displacement within four percent for tho Coast Guard designs 
tried. 
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CHAPTER I 



IHTRODUCTIOK 

Although digital conputers have been used by engineers 
for more than twenty years, their potential as design 
aids is only beginning to be realized. Naval architects 
made use of computers soon after they became available 
for solving specific algorithms such as displacement sheet 
calculations and speed-power estimates. Such algorithms 
are beneficial in the preliminary and later design stages. 
Hov/ever, with the growth of systems engineering in ship 
design, more emphasis is being placed on the conceptual or 
earliest phase of a design. At this stage of the design, 
the computer can be used more effectively than at any other 
stage . 

Because of the complexity of the naval architectural 
calculations involved in later design stages, the principle 
characteristics of a ship must be chosen early in the design. 
The output of the conceptual design phase is this set 
of characteristics. Without the use of a computer, the 
best that can be done is to develop a feasible set of 
characteristics when a optimum set is desired. Systems 
engineering is placing more emphasis on obtaining this 
optimum. 

To choose an optimum or near optimum set of character- 
istics, a great many feasibility studies are required. It 
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is impossible to do more than a few feasibility studies 
without the aid of computerized models because of the 
time required. At this early stage, the ship must be 
designed as a whole. This precludes dividing the work 
into a number of parallel paths with several people 
working simultaneously on a single feasibility study, 
jdach study is essentially a one man job. 

V/ith hand calculations a good engineer can perform 
about one feasibility study per day. As a result, great 
reliance must be placed on the designers intuition. In 
tho past this has led to series of designs which may or 
may not have been optimum. £ach new design was a slight 
perturbation of a previous successful design. While this 
procedure produces satisfactory designs, bad features 
tend zo be duplicated and few nev/ features are added. 

The engineer v/as not really at fault because no tool 
was available to perform all the calculations necessary 
to obtain an optimum design. 

Computers permit the mechanization of these hand 
calculations. The software required will be referred to 
here as a ship synthesis model. The model aids the designer 
by freeing him from the routine calculations which take 
up most of his time during the early phases of design. 

The algorithms referred to previously also free the 
designer from tedious calculations, 

Vi'hat makes a ship synthesis model different is that 
it goes beyond merely reducing the time and effort required 
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to do calculations previously done by hand. Using the 
model, hundreds or even thousands of feasibility studies 
can be made, iVhile this still does not guarantee that 
the optimum solution is obtained, the engineer will have 
much more confidence in his design than he would if only 
a few feasibility studies had been made. The designer 
is free to concentrate on the more important aspects of 
the early stages of ship design such as systems level 
tradeoff decisions. These include hull form, speed, 
endurance and payload options among others. These aspects 
were seldom considered in detail before the advent of 
ship synthesis models, 

Equally important is the fact that computers permit 
a consistency among feasibility studies. With several 
people doing independent feasibility studies or even 
with a single person doing them all, it is difficult to 
insure that the same assumptions will be made for each 
study and that no bias is applied. The model guarantees 
coriSisteucy, And consistency is important because without 
it tv/o feasibility studies can not be compared and an 
optimum solution is difficult to select. 

Several models have been available for a number of 

years and have proven their value. The U,S, Navy's 

(1 )* 

destroyer model, DDO?' ' , on which the program developed 
in this thesis is patterned, v/as started in the 1960*s 

^ Numbers in parenthesis refer to References at the 



end of this thesis 
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and updating work has been continuing since then. This 
was the first true "systems" model. Models also 
have been developed which attempt to model all surface 
ship types with a single program. For accurate results, 
each different ship type really requires its own model. 

In addition to the destroyer model, the U.S. Navy 
has developed individual synthesis models for aircraft 
carriers, underway replenishment ships, submarines, 
hydrofoils, and several other ship types. The U.S. Coast 
Guard presently has no model which gives good results for 
cutters. The destroyer model, DDO7, does not apply to 
ships as short as 150 feet in length which is one Coast 
Guard requirement. Some all inclusive models have been 
tried but with poor success, with predictions of displace- 
ment as much as 30 percent off when data on existing 
designs is inputed. 

The program, which will be referred to as the cutter 
model, developed in this thesis is an attempt to satisfy 
the Coast Guard's need. The cutter model applies to 
search and rescue and patrol type cutters with lengths 
between 150 and if -00 feet. It is not applicable to buoy 
tenders, tugs or ice breakers. An effort has been started 
by the Coast Guard to develop a model for buoy tenders 
and tugs but this v/ork is not yet complete. 

£ven in the early stages of design where the synthesis 
models are of use, a great many calculations are required 
to prove the feasibility of a ship design. This has led 
designers to employ a type of regression analysis where 
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data frorii existing ships is scaled up or down in order 
to estimate values for a new design. This estimation 
procedure is employed in ship synthesis models also, 

I'.ore detailed calculations could he made but at the cost 
of greatly increased computer time since a large number 
of feasibility studies will be run. 

The shortcut estimating techniques employed place 
more emphasis on the relative accuracy among studies 
rather than on the absolute accuracy of the result. 

Because of this, all estimating relationships in the 
model must shov; the dependence of the result on all 
significant input variables. At present, this is an im- 
possible task because not enough is known about this 
dependence. As models become more widely used much 
development v/ork v;ill have to be done in this area. An 
estimate of the variation with the most significant 
input variable is the best that can be done at present. 

Because the estimating relationships are rather 
crude, estimates are only made at the ship system level. 
Ship system level information consists of such elements 
as total outfit weight, total volume required for machinery 
and required electrical load for ventilation. Although 
total arrangements area is computed, no attempt is made 
to perform an arrangement of spaces or to compute trim. 

The engineer usually has sufficient flexibility in these 
areas to arrive at a feasible design if sufficient space 
is available. 
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One problem in developing a model for Coast Guard 
cutters is that there are few designs available from 
which to obtain data. This causes difficulty in deter- 
mining rational scaling laws. In many cases the use of 
naval architectural principles indicates the correct 
scaling procedure. In other cases all that can be done 
is to search for the most logical correlation. 

The design produced by a synthesis model should 
be tailored to meet the conditions required of a feasible 
solution. First there must be a balance between weight 
and displacement. Second, internal space available must 
equal internal space required, .This requirement is becoming 
very important as most Coast Guard vessels are volume 
rather than weight limited. Third, the energy available 
must equal the energy required. This includes both 
propulsion and electrical power requirements. Finally, 
the distribution of weight and volume must be such as to 
satisfy design criteria for transverse stability, girder 
strength and sea keeping. There are also other require- 
ments implicit in the solution such as hull type, material 
of construction and space requirements. 

In general, the inputs required for a ship synthesis 
model axe a description of the “payload" to be carried, 
a statement of vehicle mobility perf ordnance and a selection 
of options to be used. The payload description must be 
more detailed than a simple weight breakdown. The relative 
vertical centers of gravity, area requirements, electrical 
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pov/er requirements and cost must also be given. Payload 
for the cutter model also includes people, and thus the 
si?.e of the crew is an input. The description of vehicle 
mobility performance includes both maximum calm water 
speed and endurance speed and range requirements. Options 
include length, Cp, CxT as well as machinery type and 
margins. There is also an option in the cutter model 
to input the maximum sustained horsepower, in which case 
the program v/ill calculate the maximum sustained speed, 

A full description of the program inputs for the cutter 
moflel is contained in Chapter III (Main Program Description) 
and in Appendix A, 

The output generated by the model must be sufficient 
to demonstrate the feasibility of the design, V/herever 
possible, the output should permit the designer to verify 
the feasibility vdthout doing additional calculations. 

On a few of the studies run and in particular on the design 
selected, the design team will have to lay out a rough 
arrangements drawing to verify that sufficient space is 
'available. Sufficient information should be given in the 
output to permit this to be done easily. 

The output generated should also allow the designer 
to rapidly determine the design that best satisfies his 
optimization criterion. This criterion usually takes the 
form of a total acquisition cost or life cycle cost, 

*The definition of coefficients and dimensions used 
in this thesis may be found in Reference (3). 
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Predictions of acquisition cost are often included in a 
ship synthesis model. 

finally , the output should provide a basis for beginning 
the preliminary design phase. Generally the preliminary 
phase is handled by a number of design teams working in 
parallel. The output produced must provide information 
to initialize each of rhese parallel developments. 

Sample output listings for the cutter model are 
given in Appendix 3. 

It is important to emphasize that the cutter model 
produces a feasible solution only and not an optimum 
solution. An optimum solution implies that the best of 
a very large number of feasible ships has been chosen. 

It is felt that the design process works more efficiently 
if the designer acts as a guide for the computer rather 
than having the computer do the entire optimization on a 
siiigle run. In this way, many unnecessary designs can 
avoided in areas where designs are less attractive than 
others already tried. Such a process takes more real 
time but less computer time. It also gives the designer 
a much better feel for the sensitivity of the optimum 
to changes in characteristics. 

The remainder of this thesis is arranged as follows: 

Chapter II gives an overview of the cutter model 
and explains the major options. 

Chapter III goes into detail on each of the sub- 
routines in the model and explains the derivation of all 
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of tho empirical formulas used. This chapter may be 
skipped by anyone not interested in the details of the 
urogram. 

Chapter IV gives an evaluation of the model including 
test runs made using existing ships. 

Chapter V contains recommendations and conclusions. 

Appendices A and 3 contain input and output formats 
respectively. 

Appendix C is a listing of the program. 

The remaining appendices give background information 
used in this work. 

It is recommended that at least Chapters I, II, IV, V 
and Appendix A be read before attempting to use the program. 
Chapter III is aimed primarily at the person who wishes 
to modify the program, Hov/ever, the important aspects 
of each subroutine are covered in the introduction to 
each of the subsections of Chapter III. 
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CHAPTER II 



GENERAL PROGRAT.i DESCRIPTION 

This chapter gives an overall description of the 
solution method used by the cutter model but does not 
discuss the detailed calculations. The detailed calculations 
are explained in Chapter III, This chapter should also 
serve as a guide to understanding how each subroutine 
discussed in Chapter III fits into the total program. 

The organization of the computer program can be 
viev/ed in two different ways. The first is the program 
control organization which is shown in Figure 1, Sub- 
routine XECUTE controls the entire solution of each 
feasibility study case. Control is returned to the main 
program only v/hen data for a new case is required. 

Subroutine XECUTE calls other subroutines as needed and 
gradually fills all the data storage locations needed 
for output. This subroutine also controls the iteration 
loops for v/eight and vertical center of gravity. 

The logical organization of the program is shown 
schematically in Figure 2, This is an easier organization 
scheme to follow and shov;s more clearly how the program 
solves each input case. It will be noted that subroutine 
XECUTS has not been shown in Figure 2, As explained 
above, XECUTE merely guides the computer around the 
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PROGRAM CONTROL ORGANIZATION 



Figure 1 
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Figure 2 

LOGICAL PROGRAM ORGANIZATION 
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logical organization chart. Error messages are also 
printed by subroutine XSCUTS if any constraint of the 
model is violated. 

The program begins with data input in the MAIN 
program. The primary purpose of the MAIN program is to 
read data cards and all input is done in this routine. 

Some bookkeeping is also done in MAIN such as summing 
payload weights and costs. The final purpose of this 
routine is to estimate initial values of vertical center 
of gravity, KGi displacement, and cubic number. 

Once all data for one feasibility study has been 
read by MAIN and the initial values of KG, displacement and 
cubic number have been determined, control is transferred 
to subroutine XSCUTE. This subroutine directs the execution 
of the program until a solution is obtained or until • 
a constraint is violated. When either of these occur 
control reverts to the MAIjN program where data for the 
next feasibility study is read. If there is no more 
data the program stops. 

The first subroutine called by XECUTE is an under- 
water dimensions subroutine, UWDIM. Subroutine U^TDIM 
uses the input values of length, free surface correction, 

C^, and Gl.ySeam as well as the estimated KG and dis- 
placement to calculate the required beam, draft and 
v/aterplane coefficient, 0^,^. To the level of accuracy 
of this program these completely determine the under- 
water shape of the hull. 
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Two options are available with regard to speed and 
pov/er calculations. Under the first option the maximum 
sustained speed of the vessel is input. If this option 
is specified, XJICUTD calls subroutine HPCAIC which, 
together with subroutine CRVAL, estimates the horsepower 
required to maintain the maximum sustained speed. 

Option two allows the user to specify the maximum 
sustained horsepower. In this case the program will 
compute the maximum sustained speed. Subroutine HPCALC 
has been written to take speed as an input and return 
a value of horsepower. Therefore, an iterative procedure 
is required if this second option is used. Control of 
this iteration takes place in subroutine XECUTE. 

Regardless of which of the above options is used 
the endurance speed is input. Subroutine HPCALC is 
called to determine the power required at endurance speed. 

A Taylor Standard Series horsepower estimate is 
used in subroutine HPCALC. Stored values of C^. are 
used. The purpose of subroutine CRVAL is to choose the 
correct C^, value from a stored array. Reference (7) 
explains this method for estimating horsepower. 

The electrical power routine, SPLAMT, is the next 
one called. This routine determines the maximum electrical 
load at sea, the average electrical load and sizes the 
generators. The average load is needed for estimating 
generator fuel consumption. The input electrical loads 
along with crew size and cubic number are used in this 
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routine . 

The weight of fuel, lubricating oil and potable 
water is next calculated in subroutine MACHLQ. The 
endurance raiige given as input together with crew size 
and the endura^ice power calculated earlier are used by 
this routine. Lub oil requirements are either input, 
in the case where horsepower is specified, or are estimated 
based on specified machinery type. Potable v/ater require- 
ments depend entirely on the crew size. If make-up feed 
water is required it should be specified as a cargo load 
item. The specific fuel consumption is either input, in 
the case where horsepower is specified, or is estimated 
based on machinery type and horsepower. 

Subroutine M33IZ3 is next called if the sustained 
speed is specified. This routine estimates the machinery 
box length and minimum depth. If the horsepower input 
option is specified these dimensions are input. 

The next subroutine is one of the most important 
in the program. Subroutine YOLUIviE determines both the 
sheer line and the deckhouse volume. These complete 
the definition of internal volume and arrangements area 
since the underwater volume has been determined in sub- 
routine UTOIM, Because most ships for which this program 
applies are volume rather than weight limited, it is 
important to insure that there is sufficient space for 
all ship functions and equipment which must be installed. 

The volume is divided into three categories. First, 
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there is the engine room. The dimensions for the engine 
room have been estimated earlier. Second is tankage 
volume and finally arrangements volume. Arrangements 
volume is more useful if converted to an arrangements 
area. 

As a first step in determining the space available 
in each of these categories, the sheer line must be deter- 
mined, Subroutine V0LUi-i2 calls subroutine SHEER to 
determine an acceptable sheer line. The sheer line must 
meet several criteria including adequate freeboard forward 
and aft, adequate depth- amidships for the engine room 
and for adequate strength, and an acceptable sheer line 
curvature , 

V/ith the sheer line fixed, the total volume of the 
hull is estimated. The volum.e of the engine room is 
also estimated and then subtracted from the total volume. 
Next, an estimate is made of the fraction of the remaining 
volume which can be used only for tankage because it is 
unsuitable for use as arrangements volume, 

A required tankage volume is then estimated based on 
the weig'its of fuel, lub oil, potable water and aircraft 
fuel. This required tankage volume is compared to the 
volume which must be used for tankage. The larger of the 
tv;o is subtracted from the hull volume that remains. The 
remaining hull area is available for use as arrangements 
volume. This volume is converted into an arrangements area, 

A required arrangements area is computed based on 
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empirical data from past designs and on input area values. 
£y adjusting the size of the superstructure or by adding 
a raised deck, the available arrangements area can usually 
be made to equal the required arrangements area. If 
too large a deckhouse is required or if a raised deck 
longer than the length of the ship is needed, an error 
message is generated and a new input case is tried. 

Once the volume of the hull has been determined, 
the weight of all the light ship weight groups and of the 
load items can be calculated. This is done in subroutine 
V/alGHT. Although many of the items calculated by this 
routine are later printed as output, the variable of 
immediate importance is the full load weight, WFULLD. 

This value is compared with the initial guess of displace- 
ment, If there is a difference of greater than 2 tons 
betv/een the two, a new estimate of displacement is made 
and the program returns to subroutine UWDIi<i for a new 
ei.timatc of dimensions. Once a balance between WFULLD 
and displacement is achieved the program proceeds to the 
- estimation of vertical center of gravity. 

Subroutine VCG makes this estimate. Most centers 
are taken as a fraction of the average depth of the hull 
or of the amidships depth. The estimated vcg is compared 
to the initial guess and a new KG estimate is made if 
there is no agreement v/ithin a tenth of a foot. Once 
a balance is obtained, the program proceeds to a cost 
estimate . 
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The lead ship cost is estimated using a Coast Guard 
cost estimating procedure knov\Ti as "Flanagan's r.ethod"(^^ 
This procedure has been programmed as subroutine COST. 

An estimate of the speed the vessel v/ill be able 
to sustain in the North Atlantic Ocean is made in sub- 
routine oFA5?D. This routine is based on work done at 
tJie Naval ohip Engineering Center as reported in References 
(4) and (5). 

Finally, subroutine OUTPUT is called. The purpose 
of this routine is to print the values determined by the 
other routines. The input data is also printed so that 
a compact record of the feasibility study is produced. 

After the output is printed, the program returns to the 
iNAIN' program to read more data aind start a new case. 

This concludes the general description of the program. 
A much more detailed description is given in Chapter III. 
The reader v/ho is not interested in the details of the 
program can skip Chapter III with no loss of continuity. 



CHAFTiiR III 

nSTAILIiD D23CRIFTI0K OF THE PROGRAI-.: 

3.1 Introduction 

This chapter is composed of sixteen sections including 
this introduction. The other fifteen sections each 
explain in detail one subroutine of the cutter model. 

For the location of these subroutines in the overall 
program refer to Figure 2, Logical Program Organization. 

An attempt has been made to keep each section in the 
same format. A brief description of the routine will 
be given first followed by a list of required inputs. 
Statements will then be explained one by one. A list 
of the outputs produced by the subroutine will be given 
next. Finally, a nomenclature list will be given. The 
nomenclature list applies only to the subrout.ine being 
considered. In some cases the same variable will have 
different names in different subroutines. 

In both the input and output listings a code is 
employed to indicate whether the variable is transferred 
by a labeled common statement, LC(AA), or as a subroutine 
call dumray argument, 3CDA. The letters in parenthesis 
after the labeled common code indicate which labeled 
common was used to transfer the argument. 
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3.2 MAIN Program 

3.2.1 introduction 

The MAIN subroutine performs several functions. 

First and most important is the reading of all input 
data used by the model. Each feasibility study case 
starts with the reading of one or more data cards by 
the MAIN program. The actual execution of each case 
is controlled by subroutine XECUTE. This routine is 
discussed in Section 3 * 3 * However, once execution on a 
case is complete, control is returned to the MAIN sub- 
routine. A check is then made to see if there is more 
input data. If not the MAIN subroutine directs the 
computer to stop. 

Several bookkeeping tasks are performed in this 
subroutine also. These include such jobs as adding 
together the deck area and electrical load requirements 
of the individual payload items. The final function 
performed is the setting of initial values for displace- 
ment, KG and cubic number. These values are used for the 
. first iteration and are then replaced by better estimates. 

3.2.2 Inputs 

Block Data subroutines have not been used in this 
model. As a result there are no variables with assigned 
values at the beginning of the MAIN program. All of the 
variables required are read from data cards or assigned 
values in the routine. 

A flow chart for this routine is given in Figures 
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3a and 3t>. A listing of the program is given in Appendix 
C. The nomenclature list at the end of this section 
gives a description of each variable name used. 

3 . 2.3 Statement Descriptions 

A line by line discussion of the program statements 
follows. The REAL and COMMON statements have not been 
listed but may be found in Appendix C. 

REAL Statements - An explicit type definition has 
been employed so that the variable names will be more 
descriptive of their actual meanings. 

C0I/1I40N Statements - Because a great many variables 
are required to define each case , COMMON statements have 
been used to transfer values between subroutines. This 
also saves computer storage space. All COMMON statements 
used are labeled common. 

Cr array input; 

READ(«,120} CRl 

120 F0RiMAT(16F3.0) 

KEAD(8,121) CR2 

121 F0RI4AT(9F5.0) 

READ(b,l22) CR3 

122 F0RMAT(6F3.0) 

The three arrays which contain the Taylor’s Standard 
Series Cj, values are read only once at the beginning of 
the program. 

INDEX input; 

2 READ( 8,124) INDEX 

124 FORMAT(Il) 

The variable named INDEX has two functions. First, 
as described below, it tells the program when to stop. 
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Figure 3b 
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It is also used to discern between full hull ships and 
fine hull ships for the initial guess of displacement. 

Stop execution: 

IF(INDEX.EQ.O) STOP 

This is the only place where the program is directed 
to stop execution. However, if a data card is left out, 
the program will stop when there is no more data cards 
to read. 

Heading input: 

READ( a , 123 ) ( HEADd ; , 1=1 , 20 ) 

123 F0RMAT(20A4) 

Variable HEAD contains any alphanumeric information. 

It will be printed without change at the beginning of the 
output for all test cases using the data which follows, 

NO ARM input: 

R£AD(a,luOJ NOARM 

100 F0RMAT(I2) 

NOARM is the total number of armament, aircraft 
and cargo payload data ceirds which follow, A maximum 
value of 20 is allowed. 

Armament, aircraft and cargo input: 

READ(«,101)((TITLE(J,I),J=1,5).CARG0W(I),CARG0C(I), 
GR7WT( I ) , GR7CG( I ) , AInOWT( I ) , AI^0CG(I ) , ACWT( I ) , 
ACCG(I) ,ELLD(I) ,DHAR(I),HLAR(i; ,GR?CST(I),I= 

1 , NOARM) 

101 F0RMAT(5A4,12F5.0) 

Up to 20 armament, aircraft and cargo payload descrip- 
tion cards are allowed. The data may be lumped together 
and appear on only one card or it may be broken up in any 
way desired. However, if a weight is given, a vertical 
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center of gravity must also be given on the same card 
immediately following the weight. The first 20 spaces 
on each card are reserved for an alphanumeric description 
of the payload items on the card. This is used as a line 
title in the output listing. Of the four weight categories 
on the card only the group 7 weights are included in the 
light ship weight. From a calculation standpoint it is 
unimportant whether ammunition or aircraft are listed 
separately or considered cargo. However, if they are 
listed as cargo they will be classified incorrectly in the 
output listing. 

All vertical center of gravity values should be 
inputed as a multiplier of the hull depth at amidships. 

This number is easier to estimate than the actual center 
in feet. As an example, if a payload item has its center 
10 feet above the baseline and the depth at amidships 
is 20 feet the value entered is 0.5* 

The deckhouse area should include only the area 
that is required to be in the deckhouse. If the area 
could be in either the deckhouse or the hull it should 
be included in the hull area. The group 7 cost should 
include only the cost of materials and should be estimated 
for the time when the ship will be built. Because of 
format limitations the value input should be the true 
cost divided by 1000. All weights are in tons, all areas 
in square feet, all electrical loads in KW and all costs 
in dollars. This applies to all the input items. 
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Summing the payload data: 

WCARG0=0. 

WTG7=0. 

WTAT<M0=0. 

WTAC=0. 

EL0AD=0. 

BLOAD=0. 

AREADH=0. 

AREAHL=0. 

CSTG/=0. 

DO 10 1=1, NOARM 
WCARG0=WCARG0+CARG0W( I ) 

WTG7=WTG7+GR?WT( I ) 

WTAMM0=WTAMI«:0+;uM0WTC I ; 

WTAC=WTAC+ACWTCi; 

BL0AD=BL0AI>4-ELLDa ; 

GR7csT(i ;=GR7CST(i ;nooo. 

CSTG7=CSTG7+GR7CSTU ) 

AREADH=AREADH+DHAR ( I ) 

10 AREAHL=AREAHL+HLAR(i; 

Only the totals or the payload items entered are 
needed Tor calculations. These totals are computed in 
the above statements. 

NOELT input: 

READ(«,102; NOELT 

102 F0RMAT(I2) 

This variable serves the same purpose as NOARM. A 
maximum value of 20 is allowed. 

Electronics input: 

READ(«,103)( (TITLE( J,I+20) . J=1 ,5),GR4WT(I),GR4CG(I), 
ELD(I) ,DHA(I) ,HLA(I),GR4CST(I),I=1,N0ELT; 

103 FORMAT! 3A4,6F5.0) 

The comments made about armament input also apply 
here. Group 4 weights are included in the light ship 
weight. All items normally in group 4 should be included 
since no allowance is made in the program for navigation 
equipment. 
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Summing electronics input data: 

WTG4=0. 

CSTG4=0. 

DO 11 I=1,N0ELT 
WTG4=WTG4+GR4WT( I } 

ELOAD=ELOAI>+SLD(I ) 

GR^^CSTd )=GR^CST(I )*1000. 

CSTG4=CSTG4+GR4CST( I ) 

AREADH=AKEADK+DHA( I ) 

11 AREAHL=AREAHL+HLAU; 

The comments abouc summing armament data also apply 

here. Note that the area in the hull and deckhouse is 

a combined totaj. or all the payload inputs, both electronic 

and armsiment. 

Storing initial costs: 

G4CST=CSTG^ 

G7CST=CSTG7 

Because the value of variables CSTG4 and CSTG7 are 
modified during the execution of the program, a record 
is required of their initial values so that they will 
be available if another test case using the same payload 
is run. 

Manning input: 

READ(b,104; NOOFF.NOCPO, NOCREW 
104 FORMAT! 313) 

N0FF=FL0AT(N00FF) 

NCP0=FL0AT(N0CP0) 

NENL=FL0AT( NOCREW) 

The number of officers, chief petty officers and 
crew members are input and then converted to real values 
so that they can be used in non-integer equations. 

Cost indices input: 

READ! d , 1 i 0 ) DOLHR , GlIND , G2IND , G3IND , G3IND , G6IND 
110 FORMAT! 6F5.0) 
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The labor rate for the year the ship will be built 
must be estimated as must the inflation corrections for 
material costs. These corrections are for the period 
from 19^9 to the year the ship is to be built. Section 
3,14 gives some guidance as to the values to use. The 
indices can also include special material costs such as 
making the hull out of IiY-BO steel instead of mild steel. 

Military mission consumables input* 

READ( « , 1 05 ;ENDDAY , WACFUL ,CACFUL 

103 F0RMAT(3F5.0) 

The number of endurance days for dry provisions should 
be given. All other endurance days are taken as a suitable 
fraction of this number. The weight and center of gravity 
of aircraft fuel must also be given, 

HPMS=0, 

This is done so that there is no difficulty in 
calling subroutine XECUTE. 

JOPT input* 

READC»,106) JOPT 

106 F0RiMAT(I2) 

IFIJOPT.EQ,2)GOTO 1 

This variable denotes the option to be used in entering 
vehicle performance data. If JOPT equals 1, the maximum 
sustained speed is input. If JOPT equals 2, the maximum 
sustained horsepower and machinery characteristics are input. 

Vehicle performance input for JOPT equal 1* 

RE AD( a , 1 07 ) VSUS .VEND , RGEND , PCE , PCM 

107 F0RMAT(2F3.0,F10.0,2F5.0} 

GOTO 9 
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The maximum sustained speed is input along with the 
endurance speed and endurance range. The propulsive 
coefficients at the endurance speed and at maximum speed 
must also be given. 

Vehicle performance input for JOPT equal Zi 

1 KE AD( « , 1 0 b ) HPI.IS , DHI.I , XLM , WTG2;SFCMHP , SFC HHP , VEND , RGEND , 
PAL0AD,PCE,PCM,CSTG2,WTL0 
10«FOKMAT(7F5.0,F10,0,5F5.0) 

HPr-is=HPr^iE*iuoo, 

CSTG2=CSTG2*1UU0, 

GR2CST=CSTG2 

This option is used primarily when a specific machinery 
plant is desired. Besides the maximum sustained horsepower, 
the length and minimum depth of the engine room must be 
given. The total group 2 weight and specific fuel con- 
sumptions at maximum and half power are required. Tne 
endurance speed and range must be given as before. Pro- 
pulsion auxiliaries electrical load and propulsive coef- 
ficients at full and endurance speed are required as are 
the group 2 material cost and lub oil weight. 

To fit into the format used, the horsepower and 
cost values must be divided by lOUU, The cost of group 
2 must be stored for later use as were the costs of group 
4 and group 

Options input: 

y RE AD ( b , 1 09 YPE , FRSC , LEN , C P , C X , GmUL , DBI.1AR , DELCF 

109 FORrvIAT(I5.6F5.0,FlU,o; 

Machinery type must always be specified even if the 

option of inputing machinery characteristics was used. 

If machinery information was given, any number except 
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zero can be used for MTYPE. A value or zero will cause 
the program to assume the next input value is INDEX. 

This can be used to return and stop the program or to 
input a new set or data lor another ship. If the machinery 
particulars have not been specified, MTYPE equals zero 
has the same effect. However, MTYPE equals 1 will cause 
the program to assume a diesel plant is to be installed. 

A value of MTYPE equsil to 2 will result in calculations 
I or a CODOG plant. These are the only two machinery 
types for v/hich estimating relationships have been included 
in the program, if MTYPE is any other number, the results 
are unpredictable, 

A free surface correction can be specified. This 
value is interpreted as a rise in the full load KG of the 
ship in an amount equal to FRSC in feet. A negative value 
can also be used. In this case FRSC can be interpreted 
as a lowering of the average KG value built into the 
program. 

The length between perpendiculars is given next. The 

values for Cp and C^ have the standard naval architectural 
( 3 J 

definitions. GMT/iUL if the fraction GM/Heam, A value 
of 0,1 is commonly used. The next variable, DBMAK, is 
a design and builders margin. This is given as a fraction, 
i,e., a 10 percent margin is input as 0.1. DELCF is the 
ACj- value used in the calculation of horsepower. A value 
of 0.0004 is common. 

IF(MTYPE.EQ.o; GOTO 2 



One or more option input cards may follow the first 
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The program will assume that each one is a new input 
case with all other data the same. To input an entirely 
new set of data, an option input card with a value of 
MTYPS equal to zero should be included. No other values 
need be given on the card. This will cause the next 
input card to be interpreted as a new value for INDEX. 

DBKAR=1.+DBMAR 

CSTG2=GR2CST 

CSTG4=G4CST 

CSTG7=G7CST 

DBXAR is changed to a multiplier for weight and the 
input costs are reset to their initial values. 

Initial KG estimate: 

KGTRY=LEN/( .046^*LEN+4.»25; 

KGTKY=KGTRY-2. 

The initial KG value is estimated based on the values 
for current Coast Guard vessels. These are shown in 
Figure 4. Reducing the initial estimate by 2 feet was 
found to be advantageous. It helps to insure that the 
first iteration will not exceed the KG constraints of 
the program. 

Initial displacement estimate : 

DPTRY=( .006562-. UOllbfci*LEN/l 00. )»LEN**3/35. 
IF(INDEX.EQ.2)DPTRY=.0025*LEN**3/35. 

Two initial values of displacement are available. 

For an INDEX value other than 2 the displacement is taken 
as that of conventional Coast Guard vessels. This 
displacement curve is shown on Figure 4. However, if the 




(dan)/ 4 



Figure 
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vessel has a high speed-length ratio, greater than 1,3, 
it will not be possible to find the horsepower in some 
cases if the conventional ship curve is used. . This is 
due to the limitations on the Cj. arrays. By specifying 
an INDEX equal to 2 the program assumes a fine form at 
least for the first iteration. 

Cubic number estimate: 

CN=DPTRY*160 

An initial cubic number is assumed for use in estmating 
the electrical load on the first iteration. Values for 
current Coast Guard vessels range from 140 to IdO times 
the displacement, 

CALL XECUTE(HPIvlS) 

Control is transferred to subroutine XECUTE to 
calculate the vessels characteristics. 

GOTO y 

END 

The GOTO statement causes the program to loop back 
and try to read another options input data card. The 
END statement stops compilation of the routine. 



3.2.4 List of Outputs of MAIN Subroutine 
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ACCG(20) LC(EE) 

ACWT(20) LC(EE) 

AMCCG(20) LC(EE) 

AM0WT(20) LC(EE) 

AREADH LC(6B) 

AREAHL LC(3B) 

BLOAD LC(DD) 

CACFUL LC(HH) 

CARG0C(20J LC(EE) 

CARG0W(20) LC(EE) 

CN LC(BB) 

CP LC(GG) 

CRKII 52 ) LC(TT) 

CR2(972) LC(TT) 

CR3(216) LC(TT) 

A A 

CSTG2 LC(MM) 

CSTG4 LC(MM) 

CSTG7 LC(MM) 

CX LC(GG) 

DBMAR LC(WW) 

DELCF LC(CC) 

DHA(20) LC(KK) 

DHAR(20) LC(KK) 

DHm!! LC(JJ) 

DOLHR LC(NN) 

DPTRY LC(GG) 



ELD(20) LC(KK) 

ELLD(20) LC(KK) 

BLOAD LC(DD) 

ENDDAY LC(AA) 

FRSC LC(G6) 

GlIND LC(NN) 

G2IND LC(NN) 

G3IND LC(NN) 

G5IND LC(NN) 

G 6 IND LC(NN) 

GMMUL LC(GG) 

GR4CG(20) LC(BE) 

GR4CST(20) LC(KK) 

GR4WT(20) LC(EE) 

GR7CG(20) LC(EE) 

GR7CST(20) LC(KK) 

GR7WT(20) LC(EB) 

HEAD(20) LC(FF) 

HLA(20) LC(KK) 

HLAR(20) LC(KK) 

HPMS SCDA 

JOPT LC(CC) 

KGTRY LC(GG) 

LEN LC(BB) 

MTYPE LC(CC) 

NCPO LC(AA) 



List of Outputs of MAIN Subroutine (cont) 






NENL LC(AA) 

NOARM LC(FF) 

NOELT LC(FF) 

NOFF LC(AA) 

PALOAdT? LC(DD) 

PCE LC(CC) 

PCM LC(CC) 

RGEND LC(CC) 

SFCHHpTT LC(CC) 

SFCMHP LC(CC) 

TITLE(5t^0) LC(FF) 



VEND 




vsust 




WACFUL 




WCARGO 




WTAC 




WTAMMO 




WTG2t? 




WTG4 




WTG7 




WTLO 




XLM 





* Output only if JOPT * 1 
♦♦ Output only if JOPT * 2 



I 

i 
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.2.5 Nomenclature List 



ACCG(20) 


aircraft vcg/DlO, DIO » depth amidships 


ACWT(20) 


aircraft weight, tons 


AM0CG(20) 


ammunition vcg/DIO 


AM0WT(20) 


ammunition weight, tons 


AREADH 


total input deckhouse arrangements area, sq ft 


AREAHL 


total input hull arrangements area, sq ft 


AVGKW 


average electrical load, KW 


B 


maximum beam, feet 


BLOAD 


total armament electrical load, KW 


CACFUL 


aircraft fuel vcg/DIO 


CARG0C(20) 


cargo vcg/DIO 


CARG0W(20) 


cargo weight, tons 


CCARGO 


total vcg of cargo inputs, feet 


COST 


construction services cost, dollars 


CGAC 


total vcg of aircraft, feet 


CGAMMO 


total vcg of ammunition, feet 


CGGREW 


total vcg of crew, feet 


CGFUEL 


total vcg of ships fuel, feet 


CGLO 


total vcg of lub oil, feet 


CGPE 


total vcg of personal effects, feet 


CGPS 


total vcg of personnel stores, feet 


CN 


cubic number, cubic ft 


CP 


prismatic coefficient 


CRKII52) 


1 Cp arrays for Taylor's Standard Series 


CR2(972) I 

CR3(216) ^ 


, resistance estimate 



I 



Nomenclature List (cont) 
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CSTGl 


total cost of weight group 1, dollars 


CSTG2 


total cost of weight group 2, dollars 


CSTG3 


total cost of weight group 3# dollars 


CSTG4 


total cost of weight group 4, dollars 


CSTG5 


total cost of weight group 5» dollars 


CSTG6 


total ccjt of weight group 6, dollars 


CSTG7 


total cost of weight group 7, dollars 


CWP 


waterplane coefficient 


CX 


midship section coefficient 


DAVG 


average hull depth, feet 


DBMAR 


design and builders margin 


DOST 


design cost, dollars 


DELCF 


used in resistance calculation 


DHA(20) 


electronics deckhouse area, sq ft 


DHAR(20) 


armament deckhouse area, sq ft 


DHM 


minimum machinery box depth, feet 


DHV 


deckhouse volume , cu ft 


DOLHR 


labor rate, dollars/hour 


DPTRY 


displacement estimate, tons 


ELD(20) 


electronics electrical load, KW 


ELKW 


generator rated load, KV/ 


ELLD(20) 


armament electrical load, KW 


ELOAD 


total electronics electrical load, KW 


ENCVOL 


total voliune of hull and deckhouse , cu f t 


ENDDAY 


number of endurance days for dry provisions 


FRSC 


free surface correction, feet 



Nomenclature List (cont) 
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GlIND 

G2IND 

G3IND 

G4CST 

G5IND 

G6IND 

G7CST 

GMI4UL 

GR2CST 

GR4CG(20) 

GR4CST(20) 

GR4WT(20) 

GR?CG(20) 

GR7CST(20) 

GR7WT(20) 

H 

HEAD(20) 

HLA(20) 

HLAR(20) 

HPI-TS 

INDEX 

JOPT 

KGTRY 

LEN 

LRD 



cost index for weight group 1 

cost index for weight group 2 

cost index for weight group 3 

total material cost weight group 4, dollars 

cost index for weight group 5 

cost index for weight group 6 

total material cost weight group 7» dollars 

G^V^Beam ratio 

total material cost weight group 2, dollars 
group 4 vcg/DlO 

material cost weight group 4, dollars 
group 4 w*eight, tons 
group 7 vcg/DlO 

material cost weight group 7» dollars 
group 7 weight, tons 
draft, feet 

descriptive heading (alphanumeric) 
electronics hull area, sq ft 
armament hull area, sq ft 
maucimum sustained shaift horsepower 
control variable 

input option variable for machinery 
initiail KG estimate, feet 
length between perpendiculars, feet 
raised deck length, feet 
miscellaneous costs, dollars 



MOST 



1 









Nomenclature List (cont) 
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MTYPE 

NCPO 

NENL 

NOARM 

NOELT 

NOFF 

NOCPO 

NOCREW 

NOOFP 

PALOAD 

PCE 

PCM 

RGEND 

SFCHHP 

SFCMHP 



SHPE 

SHPM 

TITLE(5.40) 

TOTCST 

VEND 

VSUS 

WACFUL 

WCARGO 

WFULLD 



machinery type option 

number of chief petty officers (real) 

number of enlisted men (real) 

number of armament input cards 

number of electronics input ceirds 

number of officers (real) 

number of chief petty officers (integer) 

number of enlisted men (integer) 

number of officers (integer) 

propulsion auxiliary electrical load, KW 

propulsive coefficient at endurance speed 

propulsive coefficient at maximum speed 

endursince range , nautical miles 

specific fuel consumption at half horsepower, 

Ibs/SHP-hr 

specific fuel consumption at maximum 
horsepower, Ibs/SHP-hr 
shaft horsepower at endurance speed 
shaft horsepower at maximum speed 
line titles . (alphanumeric) 
total lead ship cost, dollars 
end\irance speed, knots 
maiximum sustained speed, knots 
weight of aircraft fuel, tons 
total cargo weight, tons 
full load weight, tons 



Nomenclature List (cont) 



WLSHIP 

WTAC 

WTAMMO 

WTCREW 

WTFUEL 

WTGl 

WTG2 

WTG3 

WTG4 

WTG5 

WTG6 

WTG7 

WTLO 

WTPE 

WTPS 

XL>1 



light ship weight I tons 
total aircraft weight, tons 
total ammunition weight, tons 
total crew weight, tons 
total ships fuel weight, tons 
total group 1 weight, tons 
total group 2 weight, tons 
total group 3 weight, tons 
total group 4 weight, tons 
total group 5 weight, tons 
total group 6 weight, tons 
total group 7 weight, tons 
weight of lub oil, tons 
weight of . personal effects, tons 
weight of personnel stores, tons 
machinery box length, feet 



( 
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3.3 Subroutine XECUIS 
3 . 3.1 Introduction 

Subroutine XECUTE controls the execution of each 
trial case. All of the remaining subroutines except 
two are called by this routine. Error messages are also 
printed by this routine whenever a constraint of the 
program is violated. 

In three cases, iteration control statements have 

been included in XECUTE. The first set of these statements 

controls the convergence of the speed-power estimate 

when maximum sustained horsepower is specified. A modified 

( 2 ) 

Newton-Raphson technique' is used with the slope ap- 
proximated by a secant to the speed-power cuirve. 

The second set of iteration control statements 
controls the convergence between estimated displacement 
and full load weight. A similar modified Newton-Raphson 
technique is used. 

The last set deals with the convergence of the 
vertical center of gravity estimate. For this set there 
is a tendency for the values to oscillate and a Newton- 
Raphson technique could not be used. The average of the 
estimated and calculated values of KG is used instead. 

While this eliminates the oscillations, convergence may 
be slow in some cases. 

A flow chart for this subroutine is shown in Figures 



5a » 5c » and 5 d 
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XECUTE SUBROUTINE 
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FLOW CHART 
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Figure 5a 
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Figure 5b 
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Figure 5c 
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Figure 5d 



3. 3*2 Input List for Subroutine XECUTE 
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DPTRY LC(GG) 

HPMS SC DA Only if JOPT » 2 

JOPT LC(CC) 

KGTRY LC(GG) 

LEN I.C(BB) 

MTYEE LC(CC) 

PCE LC(CC) 

PCM LC(CC) 

YE LC(CC) 

VS US LC(CC) Only if JOPT *= 1 

WFULLD LC(WW) 



3 . 3.3 Statement Descriptions 

REAL and COMMON statements will not be listed with 
individual subroutines but may be found in Appendix C. 

EXCKGaO. 

JKDEXol 

LNDSX-1 

These statements initialize values for the excess 
KG variable, the VCG loop counter and the weight loop 
counter respectively. 

7 CALL UV/DIM(R) 

This subroutine is called to calculate the beam 
and draft of the ship. The variable R transfers the 
value of any excess KG. 

IF(LEN.LT.O.) GOTO 999 

If a beeua to draft ratio of greater than 4.0 is 
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required to satisfy the stability requirements, LEN is 
set equal to -1000. in UWDIM. This signals XECUTE to 
print an error message and return control to MAIN. 

IF( J0PT.EQ.2) GOTO 1 

The next group of statements control the speed- 
power estimate. Two procedures are used depending on 
which perf orrsiance input option is specified. 

CALL HPC/XC(VSUS,SHP,PCM) 

IF(LEN.LT.O. ) GOTO 998 

GOTO 2 

These statements are used if the sustained speed 
is specified. Subroutine HPCALC is designed to return 
a value of shaft horsepower, SHP, if speed and propulsive 
coefficient are specified. In this case VSUS is the 
speed and PCM the propulsive coefficient used, 

LEN is used here also as a signal that sui error 
message is required, A negative length, LEN, will occur 
if the limits of the Cj. array are exceeded. The ranges 
of the Cj. array are i 

.48 < Cp < ,70 
2. < B/T < 4. 
for ,001< Cy <.003 
*5 < ^ 

for ,003< <.006 

•5 < V/Vl <1.3 

The GOTO 2 statement causes the routine to skip 
the speed estimate used when JOPT « 2, 
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1 VMAX - (HPMS*200,/LEN)»*. 33333 

If the maximum sustained horsepower, HPMS, was 
specified, the following statements are used for estimating 
the maximum sustained speed. The statement given above 
provides a rough first guess at the speed, 

4 CALL HPCALC( VMAX, SHP, PCM) 

IF(LEN.LT.O. ) GOTO 998 

SHP1=SHP 

VMAX1«.95*VMAX 

CALL HPCALC(VilAXl,SHP,PCM) 

IF(LEN.LT.O) GOTO 998 
SHP2aSHP 

Two points on the speed-power curve are calculated 
in the above statements, one point at VMAX and one point 
at ,95 X VMAX, The constant, ,95» was chosen arbitrarily 
to give a secant to the curve which has a slope close to 
the actual slope at VMAX. The same limits of the Cj- array 
apply here, 

DELHP*(SHP1-SHP2)/(VMAX-VMAX1 ) 

This is the approximate value of 

VMAX1=VMAX- ( SHPl - . 8*HPMS )/DELHP 
By definition the maximum sustained speed is the 
speed that can be made when the shaft horsepower is 80 
percent of the maximum sustained horsepower. This assumes 
the ship has a clean bottom and is operating in calm 
water. Therefore, the function to be solved isi 
f(V) » SHP(V) - .8 X HPMS * 0 

The Newton-Raphson method is an iterative procedure 
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in which I 
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^k+l * ‘ f*(Vjj) 

where f'(Vk) “ 

a V |Vjj 

In the present case : 

Vj^ = VMAX 
Vjc+i « VMAXl 

f(Vj^) » SHPl - .8 X HPMS 
f'(Yjj) 3 DELHP 

DELHP is only an approximation to however, 

it produces very good convergence. 

IF(ABS(VI^iAXl-VMAX).LE.V&lAXl*.005) GOTO 3 

The iteration stops when the difference between 

any two successive speed estimates is less than half a 

percent of the best estimate. The above iteration technique 

is very fast and convergence to these limits usually 

occurs in 3 to 4 iterations. 

VMAX = VMAXl 
GOTO 4 

These statements are used if the two speed estimates 
are not within the above limit. The new velocity estimate 
is used and another iteration is run. 

3 VS US - V^UX1 
SHK.1 a HPMS 
GOTO 5 

The above pair of values are saved for output and for 
use in other subroutines. 

2 SHPM = SHP/.8 

In the case where VSUS is specified, the value SHP 



o 

is equal to SKFh\ tines .8 as before. 

IF (MTYPE.EQ.l. AND. SHHd.GT. 14000.) GOTO 997 

If a diesel plant is specified, horsepowers of 
greater than 14000 are outside the range of validity 
of the model. An error message is printed if this occiirs. 

5 CALL HPCALC(VE.SHP,PCE) 

IF(LEN.LT.0. ) GOTO 998 

SHPEoSHP 

The endurance horsepower is calculated without regard 
to fouling or sea conditions. This statement is executed 
for both performance options, 

CALL EPLANT 

This subroutine is called to calculate the average 
electrical load at sea and to size the ships generators. 

CALL MACHLQ 

This subroutine is called to calculate the weight 
of fuel, lub oil and potable water required. 

IF(J0PT.EQ.2) GOTO 6 

CALL MBSIZE(B) 

Subroutine MBSIZE calculates the length and minimum 
depth of the machinery box. If JOPT * 2 these values 
are given as input. 

6 CALL VOLUME 

IF(LEN.LT.O) GOTO 996 

Subroutine VOLUME sets the sheer line amd balances 
the required airea and available area. If insufficient 
area is available, the value of LEN will be negative. 

CALL WEIGHT 

All the weights not previously estimated oi* inputed 
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are calculated in this routine. The primary value returned 

is WFULLD, the full load weight. 

IF(ABS(WFULLD-DPTRY).LT.2. ) GOTO 8 

The calculated full load weight is compared to the 

displacement estimate made earlier. If the two do not 

agree within 2 tons, a new estimate of displacement is 

made and the program rc iums to call UWDIM. If there is 

agreement, the program proceeds to calculate vertical 

center of gravity. 

IF(LNDEX.GE.2) GOTO 9 
WTl* WFULLD 
DP1«DPTRY 
DPTRY«WFULLD 
GOTO 10 

Two values of weight and two values of displacement 
are needed for the Newton-Raphson procedure explained 
below to work. On the first iteration, the above simplified 
procedure is used. 

9 FACT0Rb(DPTRY-WFULLD-DP1+WT1 )/(DPTRY-DP1) 

DP1=DPTRY 

WTlsV.TULLD 

DPTR Y=DPTR Y- ( DPTR Y-WFULLD )/F AC TOR 
In this case t 

f(DPTRY) a DPTRY - WFULLD(DPTRY) - 0 
Again a secant approximation is made to f* (DPTRY). This 
is the variable FACTOR. 

. DPTRYk - 

10 IF(LNDEX.GT.20) GOTO 995 

In case the iteration does not converge, this statement 
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provides a safety valve to stop the loop. An error 
message is printed, 

LNDEX=LNDEX+1 

GOTO 7 

A new iteration is run beginning with a call to UWDIM, 

8 CONTINUE 

LNDEXol 

CALL VCG(CFULLD) 

Once the full load weight is sufficiently close to 
displacement, the program proceeds to calculate KG. The 
weight counter is reset. The value CFULLD is the calculated 
KG. This value is compared with the previous estimate, 
KGTRY. 

IF(ABS(CFULLD-KGT8Y).LT..l) GOTO 11 

This iteration stops when the two KG values are less 
than 0,1 feot apart. 

IF{ JNDEX.GT.21) GOTO 994 

Again a safety valve is provided in case the iteration 
does not converge, 

IF ( CFULLD. LT. KGTRY) GOTO 12 
13 KGTRY=(CFULLI>+KGTRY)/2. 

JNDBX«*JND2X+1 

GOTO 7 

12 IF(R.LT..001) GOTO 13 

EXCKG=KGTRY-CFULLD 

In some cases the minimum beam to draft ratio may 
provide more stability than is required to satisfy the 
GM criterion. The ship is still feasible if this occurs, 
but may be over designed. When this occurs R will have 
a value greater than 0, In fact, it will bo equal to the 
excess KG when KGTRY is used in UV/DIH. The input value 
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of KGTRY is modified in UWDIM by adding R. Therefore, 
the value of KGTRY at this point in the program is the 
maximum value of vertical center of gravity for which 
a minimum beam to draft ratio will still be valid. If 
CFULLD is less than this value there is no point in doing 
another iteration. Instead, the program continues and 
eventually prints the excess KG available as an output. 

In all other cases a new estimate of KGTRY must be 
made. As mentioned earlier, this estimate is taken as 
the average of the initial estimate and final calculated 
value. This is done to prevent the two values from 
oscillating. As an example, in the case of the 378* 

WHEC, a KGTRY of l4 feet returned a CFULLD of l6 feet, 

A KGTRY of 16 feet returned a CFULLD of l4 feet. This 
oscillation may not occur in other cases, but the averaging 
estimate used will converge in every case, although pos- 
sibly at a slow rate, 

11 CONTINUE 

CALL COST(LEN) 

Subroutine COST calculates any remaining cost items 
not given as input. 

CALL SEASPD(AVSP) 

This subroutine estimates an average speed, AVSP, 
in the North Atlantic Ocean which includes the effects 
of sea state. 

CALL 0UTPUT(AVSP,EXCKG,2) 

Subroutine OUTPUT contains the \7rite statements for printing 
all the output and input data. The value 2 causes OUTPUT to 
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print both the input and output data. 

RETURN 

The present input case has been completed and control 
is returned to the MAIN program for more data. 

The error message statements are all similar in format 
and generally self explanatory. Only one typical one 
will be discussed. 

994 CALL 0UTPUT(0. ,0. ,1) 

Subroutine OUTPUT is called with a vsilue of 1. This 
directs OUTPUT to print only the input data. 

WRITE(5.106) 

106 FORMAT (///• NO BALA^^CE BETWEEN ASSUMED AND CALC KG 

WITHIN 0.1 FEET COULD BE I4ADE IN 20 ITERATIONS*///) 

GOTO 1000 
1000 WRITE(5.102) 

102 FORMAT (’PROGRAM PROCEEDING TO fCEXT INPUT CASE’///) 

RETURN 

END 

The program prints the appropriate error message 
and returns for more data. This concludes subroutine 
XECUTE. 



3.3*4 Output List for Subroutine XECUTE 



DPTRY 




EXCKG 


• . . . • e .SCDA 


KGTRY 




SHPE 




SHPM 
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3.3»5 Noaenclatur© List 

All variables have the same definition as given in 
the >1AIN program nomenclature list except for the following 



AVSP 


average speed in N.A. Ocean, knots 


CFULLD 


vertical center of gravity, KG, feet 


DELHP 


derivative of horsepower function 


DPI 


used to store DPTRY value, tons 


EXCKG 


excess KG, feet 


FACTOR 


derivative of weight function 


JNDSX 


vcg counter 


LNDEX 


weight coxmter 


R 


excess KG, feet 


RE 


same as RGEND in I4AIN 


SFCH 


same as SFCHHP in MAIN 


SFCM 


same as SFCI4HP in MAIN 


SHP 


dummy used for shaft horsepower 


SHPl 


dummy used for shaft horsepower 


SHP2 


dummy used for shaft horsepower 


VE 


same as VEND in MAIN 


VMAX 


dummy used for speed, knots 


VMAXl 


dummy used for speed, knots 


WTl 


used to store WFULLD, tons 



■n , LN 
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3.4 Subroutine UWDIM 
3,4,1 Introduction 

Three of the important variables determining the 
underwater hull dimensions are given as input. These are 
length, Cp and Two others, the beam and draft, are 

csdculated in UWDIM. The waterplane coefficient, Cy/p, 
is also calculated, but only as a linear function of 
Cp. 

Displacement is an input to this subroutine suid 
therefore the value of beam times draft is fixed since i 

B X T o Displacement x 35»/(Lt)p x Cp x C^) Ref, (3) 
Reference (3) niay be referred to for a detailed 
definition of the naval architectural variables used 
in this routine. 

The beam to draft ratio, B/T, is determined by the 

GM/B criterion. The GM available is given by the formula i 

GH = KB + BH - KG - Free surface correction Ref. (3) 

Values for KG, free surface correction and GH/B are 

inputs to this routine, A value for KB is estimated 

(3) 

using Morrish's approximate formulai 
KB =1 T - 1/3 X (T/2 + V/A) 
where T is the draft in feet 

V is the volume of displacement * 

B X T X L^p X Cp X , cu ft 
A is the waterplane area = B z L^p x C^, sq ft 
Substituting for V and A gives i 

KB - T X (5/6 - (Cp X 0;t)/(3 X 0«p)) 
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A value for BM is estimated using* 

BM a I*bp 7 ^ ^ 

where is a transverse moment of inertia coefficient 
which is assumed to he a linear function of C^p. 

A number of coefficients are used in this routine 
which are defined by the above equations. These are 
listed below. 

Cl = DPTRY*35./(LEN*CP*CX) « B x T 
C2 a .833 - CP*CX/(3*CWP) a KB/T 
C3 » LEN*CALPH/(DPTRY*35. ) * bm/b^ 

C4 a KG + FRSC 
C5 = GM/B a GMtlUL'^ 

Figure 7 gives a typical plot of the equation* 

KB + BM - GM a C2 X T + C3 x Cl^ - Cl x C5 

t3 T 

The curve shown was calculated using data for a 
378’ WHEC. 

The value of KG is also plotted in Figure 7» In 
the discussion below, KG refers to the value of KG + FRSC. 

A variable R is now introduced with a value of * 

R a KB + BM - GM - KG^^^ 

When R a 0, the stability criterion that GM » C5 x B is 
just satisfied. This is the desired solution. This 
point is shown as point 3 in Figure 7« 

Hov/ever, the solution is only valid provided point 3 
lies within satisfactory limits of beam to draft ratio. 
These limits are specified by the speed-power estimation 
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2, < B/T < 4. 
or since B « Cl/T 

< T < VST72 




It can be seen from Figure 7 that if the value of 
KG is lower than the value of KB + BM - GM at point 1 
there would be no draft for which R « 0. The miniauia 
GM that could be provided would still be greater than 






C5 X B. 

The solution nethod used in UWDIM proceeds as followst 

(1) The value of B/T at point 1 is calculated. 

Then a check is made to see if point 1 lies between the 
acceptable limits of B/T. 

(2) If point 1 lies within the limits of B/T, the 
value of R is calculated. If R is positive, no solution 
exists to R e 0 but the ship may still be feasible. In 
this case the minimum GM is greater than C5 x B. When 
this occurs the KG estimate is raised until R a 0 at 
point 1 amd the program returns to subroutine XECUTE. 

If it is found later that the actual KG of the ship is 
lower than the new KG estimate, the difference is printed 
as an excess KG. 

If the value of R at point 1 is negative , the program 
executes the procedure described in step 4, 

(3) If point 1 does not lie within acceptable 
limits of B/T, the value of R at point 2 is calculated. 

If R is positive at this point the procedure is the same 
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as that outlined in step 2. If R is negative, execution 
proceeds to step 4, 

(4) If the value of KG is greater than the value 
of KB + BM - GM at point 1 or point 2, the value of R 

at point 4 is calculated. If the value of R is negative 
at this point, no draft can be specified within the 
limits of B/T that will satisfy the stability criterion 
and the program will have to go on to another input case. 

If R is positive at point 4, the solution must lie 
between point 4 and point 1 or point 2 , whichever is 
lower. It is assumed that if point 1 is within the 
limits on B/T that the solution does not lie between 
point 1 and point 2 . 

(5) To solve for the draft at which R » 0, a synthetic 

division technique is used together with the Newton- 

Raphson iteration technique. This is sometimes known as 

( 2 ) 

the Birge-Vieta method. The equation for R must be 
modified slightly to use this technique. 

f(T) = T^R/C 2 = T^ - C 4 TVC 2 -C 1 C 3 TVC 2 + C 3 C 1 VC 2 * 0 
In the Birge-Vieta method three columns are used as 
shown below* 



1 

ai 

^2 



1 

b. 



1 

^2 



The values in the first column are the coefficients 
in the equation for f(T), namely* 

aj^ » -C 4 /C, - « — a _ r»_r'.3. 



' 2 * ^2 



a^ « 0, a^ * ^ 3^1 /^2 



i 



I 
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The values in coltxmns two and three are calculated 
using the equations i 
^i = + z X 

Ci « + z X 

where z is the current estimate of the draft, T. 

The new estimate of draft is them 

*new ” * * V°3 

New values are computed for columns two and three using 
Znew* procedure starts with the draft at point 4 

and continues until the difference between and z 

is sufficiently small. This iteration procedure is very 
fast and usually only three or four iterations are required* 
A flow chart for this subroutine is shown in Figures 
6a and 6b, 

3,4,2 Input List for Subroutine WDIM 



C4 




C5 




CP 




cx 




DISPL. . . . . 




FRSC 




LEN 





UWDIM SUBROUTINE 
FLOW CHART 
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Figure 6a 
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Figure 6b 
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3.4,3 Statement Descriptions 
DPTRY « DISPL/1.014 

DISPL has the value assigned to DPTRY in other 
subroutines. It is assumed that the total displacement 
is 1,4 percent greater than the molded displacement, 

DPTRY has the value of molded displacement in this sub- 
routine , 

CWPa,8065*C P+, 2645 

This relationship is plotted in Figure 8, 

C ALPH« , 1 2*0 WP- , 0414 

This relationship is plotted as Figure 8-8(b) in 
Reference (1), This relationship applies to destroyer 
ship types but was not checked against Coast Guard designs. 
C1«DPTRY*35./(LEN*CP*CX) 

C2=. 833-CP^»CX/( 3 . *CWP) 

C3»LEN*CALPH/(DPTRY*35. ) 

C4«C4+FRSC 

The derivation of these coefficients is described 
in the introduction to this section, 

BMIN=SQRT((C5+SCiRT(C5**2+12.*Cl*C2*C3))/(6,*C3)) 

HMIN=C1/BMIN 

BTR=BMIN/HMIN 

IF(BTR.LT.2..0R.BTR.GT,4, ) GOTO I 

The ratio of B/T at point 1 in figure 7 is calculated 
using these statements. The formula for BMIN is derived 
by setting dR/dB « 0 and solving for B. 
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RaC2*HMIN+C3*B!4IN**3-C4-C5*EMIN 
IF(R.LE.O) GOTO 2 
C4 sC4-FRSC+R 
B1=BMIN 
H»HMIN 
RETURN 

These statements compute the value of R if point 
1 lies within the acceptable B/T limits. If R is positive, 
the program returns the value of beam and draft at point 
1 and increases the estimate of KG by R. 

1 HMAX=SQRT( Cl/2. 005) 

BMIN=Ci/HMAX 

R=C 2 «HMAX+C 3 «BMI N**3 -C4-C 5*BMIN 

IF(R.LE.O) GOTO 2 

C4*rC4-FRSC+R 

B1=BI.1IN 

KsHMAX 

RETURN 

These values are calculated for point 2 in Figure 7* 
Actually they are slightly off of point 2 to insure an 
acceptable value in the speed-power routine. 

2 HMIN=SQRT(C1/3.995) 

BMAXaCl/HI^IN 

R«C 2 N+C 3 *BMAX**3 -C 4-C 5*EMAX 

IF(R.LT.O.) GOTO 6 

These values are calculated for point 4 in Figure 7 » 

If R is negative at this point, an error message must 

be printed. 

A(l)*l. 

A(2)=-C4/C2 

A(3)='-C1*C5/C2 

A(4)=0. 

A(5)=C3/C2*C1**3 

B(l)»l. 

C(l)»l. 

Zt.Kf.lIN 

3 DO 4 I = 2,5 
B(I)ttA(I)+Z*B(I-l) 

4 C(I)=B(I)+2*C(I-1) 




I 

I 



m 

i 



The columns of the Birge-Vieta method are calculated 
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in those equations. 

ZNEW=Z-B(5)/C(4) 

IF(ABS(ZIEW-Z).LE..01*ZNEW) GOTO 5 

Z=ZN2W 

GOTO 3 

A revised estimate of draft is calculated and the 
iteration is repeated imless the difference between the 
old and new values is less than 1 percent of the new 
estimate. 

5 H=ZNEW 

BleCl/H 

C4=C4-FRSC 

R«0. 

RETURN 

The results arc returned for point 3 in Figure 7. 

6 LEN a -1000. 

RETURN 

END 

These statements signal subroutine XECUTB to print 
an error message and proceed with a new data case. 



3.^.^ Output List for Subroutine UWDIM 



B1 LC(GG) 

CWP. LC(GG) 

H LC(GG) 



3»^.5 Nomenclature List 

All variables have the same definition as given 
in the MAIN program nomenclature except for the following i 
A(5) column 1 of the Birge-Vieta method 

B(5) column 2 of the Birge-Vieta method 




4 



'4 



Nomenclature List (cont) 



B1 


Beam* same as B in MAIN 


BMAX 


dummy value of beam 


BMIN 


diunmy value of beam 


BTR 


beam to draft ratio 


C(5) 


column 3 in of the Birge-Vieta method 


Cl 


coefficient^ see text 


C2 


coefficient, see text 


C3 


coefficient, see text 


C4 


KG or KG + FRSC, feet 




GMMUL in lUIN 


CALPH 


coefficient, see text 


DISPL 


DPTRY in MAIN 


HMAX 


dummy value of draft 


HMIN 


dummy V8due of draft 


I 


index variable 


R 


KB + BM - GM - KG - PRSC 


Z 


dummy value of draft 


ZNEW 


dummy value of draft 
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3.5 Subroutine HPCALC 
3 . 5,1 Introduction 

Subroutine HPCALC performs a Taylor's Standard 

f 7 ^ 

Series power estimation.''' The procedure used is very 
similar to the hand calculating procedure. However, two 
problems are of interest with regard to doing the calculations 
on a computer. First, the Cj, curves plotted in Reference 
( 7 ) must be converted to array form and stored in the 
program. Subroutine CRVAL, described in Section 3*6, 
chooses the correct value of Cj. from this stored array. 

Second, the surface area of the hull must be estimated. 

One approximate formula for wetted surface isi 
S = 1,7‘L*T + V/T^^^ 
where L is the ship length 
T is the draft 

T is the volume of displacement 
This form was used as a starting point and an empirical 
relationship was derived to bring the wetted surface 
more closely in line with that specified on the curves 
given in Reference ( 7 ), To do this a number of points 
were chosen on the curves given. These points were 
compared to the above formula and some correcting factors 
for B/T, Cy and Cp were added to reduce the difference. 

The formula given below gives values within 5 percent 
of the values calculated from the curves of Reference ( 7 ) 
over the entire range, 

S * (l,7‘L‘T+V/T)(,0053(B/T)^-.02(B/T)+3*Cy+,08»Cp+,926) 



This formula is probably not the best method for 
estimating the surface area. One method would be to 
convert the curves to arrays and interpolate to find the 
correct area. Such a method requires much more time and 
compute storage to produce a result only slightly better 
than this formula. Host of the error between a power 
calculated by this routine and by hand is probably due 
to error in the surface area prediction. 

A flow chart for this routine is given in Figure 9* 

3.5*2 Input List for HPCALC Subroutine 



B 




CP 




DELCF 




DPTRY 




H 




LEN 




PC 


SCDA 


VK 





3*5*3 Statement Descriptions 
DISPL=DPTRY/1 . 014 

DISPL is the molded displacement. It is assumed 

that the total displacement is 1.4 percent greater than 

the molded displacement. 

V0L.DISPL*35* 

BDR=B/H 

RN»VK*LBN*132050. 

CP=.075/(AL0G10(RN)-2. )**2+DSLCP 
SLRAT=VK/(SQRT(LEN ) ) 

CV*V0L/(LEN)**3 
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HPCALC SUBROUTINE 
FLOW CHART 
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Figure 
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These formulas calculate, respectively, volume of 

displacement, beam to draft ratio, Reynold’s number, 

Cf, speed length ratio, and C^. is the coefficient 

of frictional resistance calculated based on the ITTC line, 

IF(BDR.LT.2.,0R.BDR.GT,4.) GOTO 1000 
IF(CP.LT..48.0R.CP.GT.,7) GOTO 1000 
IP(SLRAT.LT..5.0R.SLRAT.GT.2, ) GOTO 1000 
IF(CV.LT. .OOI.OR.CV.GT.,006) GOTO 1000 
IF(CV.GT. .003.Ar^D,SLRAT.GT,l,3) GOTO 1000 

These statements check to make sure none of the 

constraints of the array are violated. If any is 

violated, the program will write an error message and 

go to the next input case, 

CALL CRVAL(BDR,CV,CP,SLRAT,CR) 

This subroutine is called to determine the value 
OR, the residuary resistance coefficient. The routine 
is described in section 3,6. 

CT=CF+CR*.001 

Sa(l ,7*LEN*H+V0L?H)*( , 0053*BDR*BI»-. 02*BDR+3. «CV+ 
,08«CP+,926) 

A«*. 00438*1. 9905’"^ 

EHPBH*A*VK**3*CT 

The total resistance coefficient, wetted surface 
• and finally the bare hull effective horsepower are cal- 
culated in these statements, 

EHPAPP= ( 1 .467- . 2*SLRAT )*£HPBH 

The effective horsepower with appendages has been 
based on data for the 378* WHEC, This data is sho\m 
in Figure 10 together with the above approximation, 
EHP».9*EHPAPP 

IF ( SLRAT . LT . 1 . 4 . AND . SLR AT . GT. . 9 )EHP« ( 1 . 74- . 6*SLRAT ) * 
EHPAPP 

IF(SLRAT.LE..9) EHP»1 ,2*EHPAPP 
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Figure 10 
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These equations approximate the EHP/BHPiaylor curv® 

for the 378 * WHEC as shown in Figure 11. 

SHP*EHP/PC 

RETURN 

Finally, the shaft horsepower is calculated using 

the input propulsive coefficient. 

1000 LEN» -1000. 

RETURN 

END 

If a constraint of the array is violated the 
value of LEN is made negative to signal subroutine XECUTE 
to write an error message. 

3»5»^ Output for Subroutine HPCALC 

The only output is SHFS, which is returned as a 
call statement dummy argument, 

3 » 5»5 Nomenclature List 

All variables have the same definition as given in 
the MAIN program nomenclature except for the followingi 
A dummy coefficient 

BDR beam to draft ratio 

CF frictional resistance coefficient 

CR residuary resistance coefficient 

CT total resistance coefficient 

CY Taylor's volxuaetric coefficient 

DISPL molded displacement 

EHP effective horsepower 

EHPAPP effective horsepower with appendages 

effective horsepower bare hull 



EHPBH 




1 
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Figure 11 



Nomenclature List (cont) 



PC 


propxilsive coefficient 


RE 


same as RGEND in MAIN 


RN 


Reynold's number 


S 


wetted surface, sq ft 


SFCH 


same as SFCHHP in MAIN 


SFCM • 


same as SFCMHP in miN 


SHPS 


shaft horsepower 


SLRAT 


speed length ratio 


VE 


same as VEND in MAIN 


VK 


speed in knots 


VOL 


volume of displacement 



\ 

i 

4 

I 

\ 

1 

4 

1 

1 

« 

I 



«> 

( 

I 

I- 

I 

I 



83 



3.6 Subroutine CRVAL 
3.6.1 Introduction 

Although this routine is physically one of the 
longest in the cutter model, the job it performs is one 
of the simplest. Its only function is to choose from 
stored data arrays the values of the residuary resistance 
coefficient, C^, which straddle the input point and 
then to interpolate to find the value of C^. at the input 
point. The input point is specified by four variables. 
The value of each of these four-bean to draft ratio, 

Cp, Cy, and speed length ratio-is specified as input 
to the routine. 

The constants stored in the three arrays which 
contain the Cj. data were determined using the curves 
given in Reference (7), Although the three arrays - 
CRl , CR2, and CR3 - are each one dimensional, the data 
they contain is four dimensional, varying with each of 
the four veiriables which specify the input point. Each 
of the arrays contains values for beam to draft ratios 
of 2,25, 3»0, and 3»75» Por each beam to draft ratio 
Cp was varied from 0,48 to O.70 in steps of 0.02. For 
each of these values of B/T and Cp, all the possible 
combinations of Cy and speed length ratio given below 
were stored. 

In array CRl is stored data for Cy values of 0.001 
and 0,002 and speed length ratios of 0,5 to 2,0 in steps 
of 0,1, Array CR2 corresponds to Cy values of 0.003, 
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0.004, and 0.005 plus speed length ratios from 0.5 to 
1.3 in steps of 0.1. Data for a Cy value of 0.006 and 
speed length ratios of 0.5 to 1.0 in steps of 0.1 is 
stored in CR3* 

In each of the three arrays, the first element 
corresponds to the lowest value of each of the four input 
variables. Speed lengcn ratio is increased most rapidly 
followed by Cp, and finally beam to draft ratio. 

The limits of the Cj. arrays given above are extended 

somewhat by extrapolation. Beam to draft ratios between 

2.0 and 4.0 are accepted. The limiting value on Cy has 
been extended by 0.001 for spaed length ratios between 

1.0 end 2,0. The maximum ranges that are accepted are 
given below. 

2.0 < B/T <4.0 
•^8 < Cp < .70 

.001 < Cy < .006 for 0.5 < V/^L < 1.3 

.001 < Cy < .003 for 0.5 < V/VE <2.0 

Sixteen values of C^. are required to straddle the 

input values of B/T, Cp, Cy, and V/^L. After these , 

sixteen numbers havo been chosen, the value of C^ at the 

r 

input point is found by interpolation, A total of fifteen 
interpolations are required, eight to select the correct 
speed length ratio, four to select the correct Cp, two 
to select the correct beam to draft ratio, and a final 
one to select the correct Cy. 

Even with the extensions listed earlier, the constraints 
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on th« combined C^, array cause the major limitations 
of the cutter model. The array data is listed in 
Appendix D. A flow chart for this routine is shown in 
Figures 12a and 12b, 

3,6,2 Input List for Subroutine CRVAL 



BTR 

CP . , 


. . . . . .SCDA 


CR1(1152)..., 




CR2(972)...., 




CR3(2l6).... 




cv 


SCDA 


VL 





3 , 6,3 Statement Descriptions 

I1(J,K,L,M)b(K-1)*384+( J-1)*192+(L-1)*16^-M 
I2(J,K,L,M)»(K-1)*324+( J-3)n08+(L-l)*9+M 
I3(J.K,L.M)«(K-1)*72+(L-1)*6+M 

These are statement functions which convert the 

four dimensions - J, K, L, and M - into a single dimension 

for use with arrays CRl, CR2, 8 uid CR3 resi^ctively, 

IF(BTR.GE,3. ) Kla 2 
IF(BTR.LT,3. ) Kl»l 
K2-K1+1 

An index is chosen in those statements ,f or the 
beam to draft ratio below the input point and above the 
input point. 

CPI *. 68 

3 IF(CP-CP1+.0001) 1.2,2 

1 CPl*CPl-.02 
GOTO 3 

2 L1=IFIX((CP1-.46)*50.1) 

L2-L1+1 
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In these statements, an array index is chosen for the 
Cp below and above the input point. The extra 0.0001 
in statement 3 and the 50.1 instead of 50 in statement 
2 are used to accovint for small errors in the computer 
calculations . 

CV1«.005 

4 IF(CV-CV1+. 00001) 5.6,6 

5 CVleCVl-.OOl 
GOTO 4 

13 J1«J1-1 
CVIbCVI-.OOI 
GOTO 14 

6 J1=IPIX(CV1*1001. ) 

IP(J1.EQ.5.AND.VL.GT.1. ) GOTO 13 
IP(Jl.EQ.2.AhD.VL.GT.1.3) GOTO 13 

14 J2=J1+1 

Arreiy indices for Cy are chosen below and above 
the input point. The statements connected with statement 
13 extend the limits of the array by extrapolation. 

VLl-1.9 

7 IP(VL-VL1+.0001) 8,9.9 

8 VL1»VL1-.1 
GOTO 7 

9 M1=IPIX( (VL1-.4)*10.1) 

M2=Ml+l 

Again, array indices are chosen, this time for 

speed length ratio values. 

J-Jl 

INDEX»1 

The value of Cy is held fixed each time through 
the following statements. When values at each of the 
two Cy values have been calculated, an interpolation 
is done to find the Cj. value at the correct Cy value, 

11 IP(J.GB.3) GOTO 10 
I»I1(J,K1.L1,M1) 

CRA»CR1(I) 
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IaIl(J,K2,L2.?42) 

CRHsCRl(l) 

GOTO 50 

Variables CRA to CRH are used for the eight Cj, values 

that remain when Cy is held constant. In this sot of 

statements array CRl is used. 

10 IF(J.EQ.6) goto 12 
I»I2(J,K1,L1.1U) 

CRA*CR2(I) 



I«I2( J,K2.L2,M2) 

CRK«CR2(I) 

GOTO 50 

This group of statements is identical to the proceeding 

group except that array CR2 is used. 

12 I=I3( J»Kl,Ll,t!l) 

CRA«CR3(I) 



I=I3( j»K2,L2,M2) 

CRH»CR3(I) 

Again these statements are identical with array 

CR3 being used. 

50 Vt.(VL-VLl)*10. 

C»(CP-CP1)*50. 

BTR1«FL0AT( (Kl-1 ) )*. 75+2.25 

B=(BTR-BTR1)/. 75 

C1»(CRB-CRA)»V+CRA 

C2*(CRD-CRC)*V+CRC 

C 3= ( CRF-CRE ) *V+CRS 

C4= ( CRH-CRG ) »V+CRG 

C5=(C2-C1)*C+C1 

C6e(C4-C3)*C-}-C3 

C7«(C6-C5)*B+C5 
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These are the interpolation statonents required 

to first calculate the value of C^, at the input speed 

length ratio, than at the input Cp and finally at the 

input beam to draft ratio, 

IP(INDEX,GE,2) GOTO 51 

A=C7 

J=J2 

INDSX=2 

GOTO 11 

When INDEX * 1 the lower value of Cy is used in the 
interpolation above. The next step is to do the interpolation 
at the upper value of C^, The value of at the lower 
Cy is saved as variable A, 

51 DsC? 

CRa(D-A)*(CV-CYl )*1000,+A 

RETURN 

END 

One final interpolation is performed to calculate 
the desired value of CR and the routine retxima control 
to subroutine HPCALC, 

3.6.4 Output for Subroutine CRVAL 

The only output is the calculated value of CR which 
is returned as a subroutine call dummy argument, 

3 . 6.5 Nomenclature List 

All variables have the same definition as given in 
the MAIN program nomenclatxire except for the followingi 
A value at lower Cy value 

B fraction of BTR interval 

BTR input value of beam to draft ratio 

lower value of bean to draft ratio 



BTRl 
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fraction of Cp interval 



Cl 


Cp at 


C2 


Cp at 


C3 


Cp at 


C4 


Cp at 


C5 


Cp at 


C6 


Cp at 


C7 


Cp at 


CPI 


lower 


CR 


Cp va; 


CRA 


Cp at 


CRB 


Cp at 


CRC 


C„ at 


CRD 


Cp at 


CRE 


Cp at 


CRF 


Cp at 


CRG 


Cp at 


CRH 


Cp at 


CV 


input 


CVl 


lower 



Cj* at correct V/7^, lower BTR, lower C 



D 

I 

INDEX 

J 

J1 

J2 

K1 



'P 

'p 



P' 

'P' 



Cp at higher V/^ and Cp» lower BTR 



ower V/Vl, higher BTR and Cp 
Lgher V/Vl» BTR, and Cp 
SLlue of Cy 
ilue of Cy 
Cp value at upper Cy value 
array index 

Indicates whether Cy is at lower or upper value 
used for either J1 or J2 
index for lower value of Cy 

index for higher value of Cy 
index for lower value of BTR 




.1 



K2 

LI 

L2 

Ml 

M2 

Y 

VL 

VLl 
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index for higher value of BTR 
index for lower value of Cp 
index for higher value of Cp 
index for lower value of V/,^ 
index for higher value of V/^ 
fraction of V/^L interval 
input speed length ratio 
lower value of speed length ratio 
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3.7 Subroutine EPLANT 
3,7.1 Introduction 

This routine porforas two functions. The first is 
to estimate an average plant electrical load for use 
in calculating fuel requirements. The second is to 
size the generators. Generator size is used primarily 
for estimating weights, but it is also a useful output 
which provides information needed for beginning the 
preliminary design phase. 

Three ships were used in developing electrical 
load relations. These were the 378* WHEC, the 210* WMEC, 
and the 150' WPB, The electrical loads for these vessels 
were categorized as shown in Table 1, This table also 
gives the loads in KW that were usod. 

The first t\ro categories, electronics load and arm- 
ament load, are inputs. Propulsion auxiliary load includes 
only equijHnent needed for the direct operation of the 
main propulsion machinery. For example, lub oil service 
pumps are included in this category but lub oil transfer 
pumps are included under other auxiliaries. Air condition- 
ing and ventilation also includes electric heaters. 

Hotel is primarily galley and lighting. 

Estimating relationships were developed for each 
of the last foxir categories in Table 1, If the propulsion 
auxiliary load is not input, a value of 1 KW is assumed 
for diesel plants and 71 KW for CODOG plants. Air 
conditioning and ventilation is calculated as a function 
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Table 1 

ELECTRICAL LOAD REQUIREMENTS 



Electrical Load 


150* WPB 


210 


• WMEC 


378* WHBC 


Electronics 


4.8 




21.99 


117.84 


Armament 


42.7 




0.0 


46.64 


Propulsion Aux. 


1.02 




.71 


71.03 


Air Cond. & Vent. 


21.1 




49.13 


308.34 


Hotel 


20.9 




35.4 


170.51 


Steering gear 


.54 


4.48 


37.38 


Aux. Machinery 


18.26 40.08 


70.72 


Shops 


1.0 


2.2 




4.98 


Other Aux. (total) 


19.8 




46.76 


113. 08 


Total 


110.32 


153.99 


827.44 



All loads are in KW 



of cubic number tines the number of accommodations. Hotel 
load is a function of the number of accommodations only. 
Other auxiliaries load is assumed to be a function of 
cubic number only. These relationships «ire plotted in 
Figures 14, 15 and l6. 

The total electrical load required is equal to the 
sun of the six categories above with a two i>ercent growth 
margin on all equipment plus a 100 percent margin on 
electronics. An additional 25 percent is included above 
the groYTth margin to allow for deterioration and starting 
loads. 
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It is assumed that three generators will be installed. 
Each generator is sized so as to be able to carry the full 



electrical load. This represents the current Coast Guard 
practice for new designs. Less capacity has been installed 
on past designs. 

The smallest generator that will be installed is 
assvuned to be 100 KW. Other possible plant sizes lie 
between 100 KW and 250 KW in steps of 50 KW, between 
250 KW and 1000 KW in steps of 250 KW, and above 1000 KW 
in steps of 500 K\7. The smallest generator that is able 
to supply the total load with the margins stated above is 
chosen. 

A flow chart for this routine is shown in Figure 1?. 
3 , 7.2 Input List for Subroutine BPLANT 



BLOAD 




CN 




ELOAD 




JOPT 




MTYPE 




NCPO 




NENL 




NOPF 




PALOAD 
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EPLANT SUBROUTINE 
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Figure 17 
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3. 7*3 Statement Descriptions 
IP( J0PT.EQ.2) GOTO 1 

The estimate of PALOAD is skipped if a value is 
given as input. 

PAL0AD=1 . 

IF(MTYPE,EQ.2) PALOAD*?!. 

Propulsion auxiliaries electrical load is calcxxlated 

for diesel and CODOG plants respectively, 

1 OALOAD*.0002695*CN 
NAC=N0FF+NCP0+rffiNL 
HL0AD=.917*NAC 

AYLOAD*15.+.376*NAC«CN/100000. 

KW»ELOAD4-BLOADfPALOAD4-OALOAD+HLOAD+AVLOAD 

AVGKW=KW-BLOAD 

KW1=1 . 25* ( 1 . 02*KW+EL0AD ) 

The individual load items are calculated and summed 



to give an estimate of the total load. The average load 
is taken as the total load minus the armament load. 



Margins are then added to the total load for sizing the 
generators. 

ELKW=100. 

2 IF(ELKW.GE.KWl) GOTO 3 
IF(BLKW.LT.2^0. ) GOT 4 
IF(ELKW.LT.900, ) GOTO 5 
ELKW-ELKW+500. 

GOTO 2 

4 ELKW=ELKW+50, 

GOTO 2 

5 ELKW*ELKW+250. 

GOTO 2 

3 RETURN 
END 



These statements choose the generator rating, ELKW, 
which is the next size larger than the total load with 
margins , 
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3,7»^ Output List for Subroutine EPLANT 



AVGKW. 




ELKW. . 





3 . 7*5 Nomenclature List 

All variables have the same definition as given 
in the MAIN program nomenclature except for the following* 



AVLOAD 


air conditioning and ventilation load, KW 


HLOAD 


hotel load, KW 


KW 


total load without margins , KW 


KWl 


total load with margins, KW 


NAC 


number of accommodations 


OALOAD 


other auxiliaries load, KW 


RE 


same as RGEND in MAIN 


SFCH 


same as SFCHHP in MAIN 


SFCM 


same as SFCMKP in MAIN 


YE 


same as VEND in MAIN 
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3.8 Subroutine I4ACHLQ 
3,8.1 Introduction 

At this point in the program, the shaft horsepower 
and size of the electricaLl plant have been determined and 
the weight of fuel, lub oil and potable water can be 
calculated. Subroutine MACHLQ performs these calculations. 
In conjunction with the fuel weight calculations, the 
specific fuel consumption, SFC, at endurance power is 
estimated. 

For diesel plants it is assumed that the SFC at 
rated power is 0.42 Ibs/shp-hr and at one half rated 
power is 0,46 Ibs/shp-hr. A linear illation is used 
for other power ratios. The SFC used for CODOG plants 
is 0.42 if the endurance power is less than 7000 SHP. 

This assumes that the diesels will be sized to be at 
maximum power at endurance speed. If the endurance 
power requirements are greater than 7000 SHP, the gas 
turbines must be used at endurance speed. Three maximum 
power levels for the gas turbines are used - 25000 HP, 

35000 HP and 50000 HP. At each of these horsepower 
levels the SFC is taken to be 0.50 Ibs/shp-hr, At half 
power, the SFC is , 595 » A linear relation is used for 
other power levels. 

If a specific machinery plant is input, the specific 
fuel consumptions at half and full power are given. A 
linear relation is used for other power levels. 
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The specific fuel consumption of the ship service 
generators is assumed to be about 0.48 Ibs/shp-hr or 
0.65 Ibs/KW-hr. The generator fuel consumption is computed 
using the average electrical load. 

Fuel consumption for hotel services is taken to be 
0.32 Ibs/man-hr. This assumes some waste heat recovery. 

The lub oil requirements for diesels is assumed to 
be 5 tons plus 1 ton per thousand horsepower. For CODOG 
plants this same relationship is used for the diesel part 
of the plant (endurance power is used in computing the 
requirements). However, additional lub oil for the gas 
turbines is included. For the gas turbines, the oil 
requirements are 3.5 tons plus 1 ton per ten thousand 
horsepower. 

Potable water requirements are taken to be 50 gal/man 
or 0,186 tons/man, A flow chart for this routine is 
shovm in Figure 18. 

3,8.2 Input List for Subroutine HACHLQ 



AVGKW 




RE 




JOPT 




SFCH 




MTYPE 




SFCM 




NCPO 




SHPE 




NENL 




SHF^ 




NOFF 




VE 





m 
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3.8,3 Statement Descriptions 
NAC=NOFF+NCPO+NENL 

The total number of accommodations is calculated. 
WTPS=.186<*NAC 

Additional items will be added to the weight of 
personnel stores later but for now it includes only the 
weight of potable water. 

IF(J0PT.EQ,2) GOTO 1 

This skips the lub oil and standard SFC estimates 
since these values eire inputs, 

V/TL0=..001*SHFK+5. 

IF(MTYPE.EQ.2) WTLOa. 0001*SHPIi+. 001*SHPE-{-8,5 

The lub oil weight for diesel and CODOG plants 

respectively are calculated here, 

IF(MTYPE.EQ,2) GOTO 2 
SFC= . 5-. 08--SHrE/SHKI 
SKPH=.SHPM/2, 

IF ( SHK4 . GT . 7000 . . AfiD . SHPE , LT . SHPH )SFC= . 5- . 1 6*SHPE/SHFM 
GOTO 3 

These statements calculate the SFC for diesel plants. 

In case the shaft horsepower is greater than 7000 it is 

assumed that four diesels will be installed. If the 

endurance power is less than half of the maximum, only 

two diesels need be in operation at endurance speed, 

2 IF(SHPE.LT.7000. ) GOTO 4 
SFC=(1.38-.38*SHPS/50000. )*,5 

IF(SHPM.LS, 25000. )SFC=^(1.38-.38^HPE/25000. )*,5 
IF(SHPM.GT. 25000., AND. SKPM.LE. 35000. ) SFC*(1,38- 
.38*SHFE/35000.)*,5 

GOTO 3 
4 SFC=,42 
GOTO 3 



For CODOG plants the SFC is calculated using the 
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statements above • 

1 SFCaSFCM+(l.-SHPE/SHI^)*2.*(SFCH-SFCM) 

If SFCH and SFCM are input, the above linear relation 

is used for calculating the SPC at endurance speed. 

3 FRMAIN=SFC«SH?E 
FRH0TL=.32*NAC 
FRGENa.65*AVGKW 
FR«FR]=1AI N+FRKOTL+FRGEN 

All of these fuel rates are in Ibs/hr. 

V/TFUEL=FR»RE/( 1788 . *VE ) 

RETURN 

END 

The factor 1?88. is derived from 

2240. *(1 ~ F .?. p ) 

Pm^^deg^^sc 

where i 

F^p * .05 is a tailpipe allowance 
Fn « 1*03 is an extra margin 
F^eg* i *05 is a degradation allowance 
Fgc « 1*10 is an allowance for sea conditions 
and fouling 

The value 2240 converts pounds to tons* 

3*8.4 Output List for Subroutine MACHLQ 

WTFUEL LC(LL) 

WTLO LC(LL) 

WTPS LC(LL) 

3*8.5 Nomenclature List 



All variables have the same definition as given in 
the MAIN program nomenclature except for the followingi 



Nomenclature 
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List (cont) 


FR 


total fuel rate in Ibs/hr 


FRGEN 


generator fuel rate in Ibs/hr 


FRHOTL 


hotel fuel rate in Ibs/hr 


FRMAIN 


main engine fuel rate in Ibs/hr 


NAC 


total number of accommodations 


RE 


same as RGfEND in MAIN 


SFC 


specific fuel consumption at endurance power 


SFCH 


same as SFCHHP in MAIN 


SFCM 


same as SFC^mP in MAIN 


SHPH 


one half maximum shaft horsepower 


VE 


same as VEND in MAIN 
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3.9 Subroutine tlBSIZE 
3.9,1 Introduction 

Very littlo data is available concerning the minimiun 
allowable depth and length of the nachinery compartment , 
As a result, only very simple calculations are made in 
this subroutine to estimate these two dimensions. 

The product of length of machinery box times the 
maximum beam of the ship was calculated for two diesel 
designs, the 210* WMEC and the 150* WP3, and for one 
CODOG plant, the 378* WHEC. A value of about 1050 sq. 
ft, was obtained for both diesel designs and a value of 
2300 sq. ft. was obtained using the CODOG plant. These 
values are divided by the beam of the ship to obtain 
the required engine room length in the program, 

A check is also made to insxire that a minimum length 
is not exceeded. This minimum length was taken to be 
28 feet or $6 feet for diesel plants below or above 
7000 BHP respectively and ^5 feet for a CODOG plant, 

A rough deck height approximation was derived using 
data from the same three vessels mentioned above. The 
formula used gives a deck height of about 1? feet for 
the 150* WPB and about 19 feet for the 210* WMEC, The 
deck height for all pleints with horsepowers greater than 
7000 is taken as 22 feet, 

A flow chart for this subroutine is shown in Figure 
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Figure 
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3.9»2 Input List for MBSIZE Subroutine 

SODA 

MTYPE LC(CC) 

SHPM LC(CC) 

3,9,3 Statement Descriptions 

IF(MTYPE.EQ,2) GOTO 1 
XU1=1050,/B 

IP(SHFil,GT.7000,,AND,XLM.LT.56. ) XLJl-56. 
IF(SKi’I/l.LS. 7000 .. AND. XUj.lt. 28. ) XUI-28. 

GOTO 2 

These statements set the engine room length for 
diesel plants. (MTYI^ « 1 for diesels) 

1 XLM=2300,/B 
IF(XU>1.LT.45.) XLM=45. 

These statements set the engine room length for 
CODOG plants, 

2 DHr.l=il. + .00157«^HFM 
IF(D}{M,GT.22. ) DKM=22. 

RETURN 

END 

The minimum deck heighli which is the same for both 
plant types, is set by these statements, 

3.9*^ Output List for Subroutine HBSIZE 

All variables have the same definition as given in 
the MAIN program nomenclature except for the following* 



RE 


same 


as 


RGEND in MAIN 


SFCH 


same 


as 


SPCHHP in MAIN 


SFCM 


same 


as 


SFCMHP in MAIN 


VE 


same 


as 


VEND in MAIN 
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3,10 Subroutine VOLUME 
3,10.1 Introduction 

Many of the ships designed today are volume rather 
than weight limited. Because of this it is important 
that a ship synthesis model determine a balance bet\«foen 
required and available volume as well as between weight 
and displacement. Subroutine VOLUME performs this function 
in the cutter model. 

The dimensions of the hull below the load waterline 
are determined in subroutine UWDIM. Only the sheer line, 
deckhouse volume and the length of raised deck, if any, 
remain to be determined to define the volume available 
in the ship. The enclosed volume of the hull and deck- 
house will bo completely determined by these variables 
since the waterplane coefficient is calculated in sub- 
routine UWDIM and since the general shape of the hull 
is implied by the estimating relationships used in the 
model. 

Deckhouse volume is constrained primarily by stability 
considerations and by external deck area roquirements. 

■ If the deckhouse becomes too high, it will not bo possible 
to make the vessel stable* Also a certain amount of 
free deck area is required in nearly all designs, restricting 
the horizontal spreading of the deckhouse. At the other 
extreme, the ship may have excessive stability if too 
small a deckhouse is installed* 

These deckhouse volvime restrictions have been included 
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in the cutter model by restricting deckhouse sise to be 
within the limits of past practice. Both em upper and 
a lower limit are placed on the size of the deckhouse 
as shown in Figure 20, The deckhouse volxxrae of Coast 
Guard vessels has been generally greater than that of 
Navy designs. Upper and lower limits of 0,00150 x and 

3 

0,001 X L , respectively, are used in the Navy’s destroyer 
model. This difference is probably due to the fact that 
the Navy requires more external deck area for armament. 

The deckhouse volumes shown in Figure 20 do not 
include uptakes. It is assumed that sufficient additional 
area is available for uptakes since uptake area was 
excluded from the data points also. However, it is 
important to check a design be sketching a layout for the 
external decks and the superstructure before acceptance. 
The volume in the hull is also estimated in VOLUME, 
The underwater volume is already available. The first 
step in determining the above water hull volume is the 
development of an acceptable sheer line. This step is 
performed by subroutine SHEER and is discussed in con- 
junction with that routine. 

Once the sheer line and with it the average freeboard 
have been estimated, the above water volume is found 
by multiplying the area of the waterplane by the average 
freeboard and by a factor to accoxmt for flare. The 

total volume in the hull is the sum of the above and below 
water volumes. 
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This total hull volume includes the volume to the 
main deck only. The hull volume cannot he reduced below 
this size but a raised deck can be added to increase the 
volume of the hull. For vessels of the size of Coast 
Guard cutters it is impractical to add more than one 
raised deck. Also the deckhouse size can be varied some- 
what to change the available volume. 

Required hull volume is broken down into three cat- 
egories. These are the machinery box* tankage volvime 
and arrangements volume. The object of subroutine VOLUME 
is to alter the available volume so that there is sufficient 
space for each of the three categories of required volume. 
The model assumes that all tankage volume and the machinery 
box volume are in the hull. The deckhouse contains only 
arrangements volume. 

The machinery box size is the first category estimated. 
This volume is found by multiplying the midships area 
by the length of the machinery box and by a machinery 
box prismatic coefficient. Machinery box volume is 
subtracted from the available hull volume to the main 
deck leaving a volume which can be used for tankage and 
arrangements. 

Some of the remaining volume cannot be xised for 
arrangements because of hull shape. The next step is to 
estimate the fraction of the remaining volume that can 
be used only for tankage. This is the minimum tankage 
volume for the ship. If loss tankage is required for 
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the liquids that must be carried, the excess volume oust 
be used for voids or cofferdams sine© it cannot be con- 
verted into machinery box space or arrangements volume. 

If more tankage is required than this minimum, some of 
the arrangements volxioe must be used for tankage. 

The weights of all the liquids which must be carried 
are calculated in subroutine HACHLQ. In VOLUME, these 
weights are converted to volumes and added together. 

A twenty percent allowance is added to this figure to 
account for required cofferdams and voids. The volume 
required for liquids is then compared to the minimiua 
tankage volume to determine the amount of arrangements 
volume, if any, that is required for tankage. 

All of the volume remaining in the hull after tankage 
volume and machinery box volume have been subtracted s can be 
used for arrangements. The entire superstructure volume 
is also available for arrangements. Arrangements space 
requirements are more commonly given as areas rather than 
8 is volumes, so the next step is to convert the available 
arrangements volume to an available arrangements area. 

For the deckhouse this is a simple task since the deck- 
house is nearly wall sided. The volume is divided by 
the deckheight of 8,5 feet to find the area. For the 
hull no such simple procedure is possible because of 
hull curvature . 

The Navy performed extensive work on this problem 
while developing their destroyer model, DD 07 , This work 
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was made use of in the cutter model not only for the 
arrangements area formulation hut also for the flare 
factor, machinery box prismatic coefficient and tankage 
volume fraction estimates discussed earlier. The relation 
ships used in the cutter model were taken directly from 
the original notes of Mr. Jim Mills who was largely 
responsible for the development work done by the Navy. 

The major deviation between the cutter model and the 
destroyer model lies in the engine room location. The 
destroyer model assumes that the engine room is centered 
at amidships but Coast Guard practice usually places the 
engine room aft of amidships. No correction is made for 
this difference and the output produced for test runs 
with existing ships indicates that there is little effect. 

By using the Navy’s formula, the arrangements area 
in the hull can be calculated. This area plus the area 
of the deckhouse must be compared to the area required 
for all the ship's functions plus the area specified as 
input. The area required for ship's functions such as 
stores area, living area and auxiliary machinery area, 
are estimated in the cutter model based on past practice. 

Data was obtained by measuring the areas of five 
past designs and assigning the areas to various categories 
This data is summarized in Table 2, Linear relationships 
were developed from this data as shown in Figures 21 
through 25. Berthing and sanitary areas are based on 
minimum Coast Guard standards. The area required for 
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CATEGORY 


378 

WHEC 


327 

WHEC 


255 

WHEC 


210 

WPEC 


150 

WPB 


Office spaces 


493 


465 


347 


130 


65 


iViessing Facilities 


2713 


2348 


1379 


1452 


339 


Crew Special 


725 


609 


717 


324 


32 


Officer Staterooms 


1157 


1338 


966 


567 


187 


Officer Sanitary 


243 


133 


120 


193 


30 


C ?0 Staterooms 


1040 


364 


336 


283 


110 


CPO Sanitary 


250 


100 


- 


75 


27 


Crew Berthing 


3641 


2655 


2568 


1658 


375 


Crew Sanitary 


716 


444 


492 


384 


90 


CO Stateroom & Cabin 


697 


511 


273 


406 


- 


CO Pantry 


58 


128 


22 


- 


- 


Commissary Stores 


600 


724 


531 


480 


35 


Other Stores 


3015 


2354 


2397 


1294 


435 


Workshops 


1161 


818 


752 


544 


- 


Passages 


2567 


2588 


2138 


739 


229 


Repair Lockers 


148 


18 


97 


110 


- 


Steering Gear < 5 : Windlass 


736 


548 


645 


270 


294 


Chain Locker 


47 


30 


•36 


42 


16 


A/C & Fan Spaces 


369 


324 


39 


165 


- 


I.C. & Gyro Room 


312 


143 


232 


180 


- 


Aux. Machinery Spaces 


I860 


1874 


4 l 4 


820 


543 


Uptakes (Hull) 
Pilothouse , Chartroom & 


462 


1216 


64 


— 


• 


C.I.C. 


1125 


856 


521 


422 


300 


Sub total 


24135 


20566 


15086 


10536 


3107 


EXTRAS 


Crew Recreation Rooms 


1445 


•• 


682 




71 


Cargo Flume Tanks 


1219 


249 


- 


- 


241 


Electronics S Stores 


l 6 l 8 


775 


1206 


319 


157 


Armament d Stores 


844 


822 


730 


280 


- 


Balloon Shelters 


512 


156 


120 


- 


- 


Mission Offices ic Labs 


746 


336 


482 


- 


- 


Mission Workshops 


- 




70 


105 




Bow Thruster 


344 


— 






• 


Total 


30663 


22926 


18376 


11242 


357^ 
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Figure 22 
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repair lockers is taken as a constant 100 square feet. 

The area for chain lockers as 40 square feet and the area 
for uptakes in the hull as either 6$ or 460 square feet 
depending on whether diesels or gas turbines, respectively, 
are installed, A constant value of nine percent of the 
total area is assuraed for passages. Passages ranged 
from 6,84 to 13*17 percent of the total area on the 
five designs used. 

All of the required areas are sunmed to detenaine 
the total required area* This area is then compared to 
the available arrangements area in the hull and super- 
structure. 

Comparisons are made in the following order. First, 
the arrangements area required to be located in the 
deckhouse is compared to the area of the largest allowable 
deckhouse. If the required area exceeds that available, 
the design is infeasible and an error message must be 
printed. Next, the total required arrangements area 
is compared to the total available area. If the required 
area is the smaller, tho deckhouse size is decreased. 

The deckhouse size used is either the size which balances 
the area requirements or the smallest allowable deckhouse 
size or the area required to be in the deckhouse, whichever 
is greatest. Should the total required area exceed 
that available a raised deck is added. 

The length of raised deck is determined by an iterative 
procedure. On each iteration the area required to be 
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included in the raised deck is coraputsd and the curve 
Bhovm in Figure 26 is used to convert the area to a 
raised deck length. This curve is based on the main 
deck of a 378' WHEC. The iteration continues xmtil two 
successive iterations fall within 5 percent of the length 
of the ship or until twenty iterations is reached. In 
the later case the design is ruled infeasible. A check 
is also made to insure that the raised deck length does- 
not exceed the length of the ship. 

If the raised deck length lies between 0,4 and 0^6 
times the length of the ship, the raised deck is extended 
to 0,6 times the length of the ship. 

A flow chart for this routine is shown in Figures 



2?a and 27b, 








3,10,2 Input List 


for Subroutine 


VOLUME 




AREADH 




LBN 




AREAHL 




LMB 




B 




MTYPE 




CP 




KCPO 




CWP 




NENL 




cx 




NOFF 




D 




WACFUL 




DMB 




WTFUEL 




H 




WTLO 




JOPT 




WTPS 
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T A R.T VOLUM B 



REAL AMD 
LASBLSD 

STA^SM£>J-rS 



KHOEX^<P. 

rJCtJL 

(f, 

Le.ou^^ 



Ax') - 

XTPAx- ( cP/VN,x,C/*) 
V/rv\8^ LA^a if/vn-jecPAi 
At^VgM 



P7= fCc>^^A) 
HV6^' L£jJn 




CAUL ^Hsee^ 

F!^ 0^ Uf,DZA 

OC£fc/rr^ 



Ojx. ((Ol^,CFOi,^Lf^2>/uS^) 
ATVcw\:x ^^4/+v'^>A\ 

V FV uJ TPC/flA>f 3. 

VL.O z tvrLo^^^, 

I/AO C 

VFCP0tLTU^9Ji^4t^,l 



U^f^tL f\/UD 
4vx^u;t v/f<P4/0«- /#2. 
/erAA 5 . 




rW>M = ^>^A + 

avaa 

ooa-ZiZM- 

LCj0^232- 



yes 







V ^TYPf3/ > 








nIo 


lA'CA; 






LY£^ 


/co/^r, \ 



irrr^ 




i 


. 

Ijo 



fOe:Kr 

PA6E 





ftTAA ■= 



cfcj ^ (o ua4joM5y 
A4_/CX>i»d 

. oolI^j44:AJ+Z3?>ieU) 
/.O^ 



/?iV/<?^-A7V0>n 
TTf AJ/dULz /4rit>^ ^ATAAi^^ 
p~. fCt>n, CP, cn,tjH«/v.e*)) 
Ai/Afts <^fiV^^>tP-/eTAA 
.o^>^^3■»«^M*Aa 

D>#AAi 0>i'J/e.g 




X.'HPDAi//,?^/V<4C +-5Ba 
C5DAi-f,Z3>»AMC 
OSADfitz. 7’I^AJOFF 
OrtDA Ts. i7>^aJ0FF 
CPo^*‘ 40,'A‘ f^ePo 
CFoAOfk Ac^O 

C^D/^--c Zi..i^AJCtJi. 

ChOA ^ S, ^fL>C/JL 



COPA r 2.37o«£:k} - 76 . 
C57D4 5 3o 0 .-►. oz/A» 
<i>i^>Az.O0'?/9^CA7 
UJ 04 




RLDA^OO, 

Cuda^4(^, 

UTdA^(P, 



Figure 27a 
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Figure 27b 
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3.10*3 statement Descriptions 
KNDEX-0 

NACaNOFF+NCPOfNENL 

deckht«o. 

LRDsO. 

LRDlaO. 

KNDEX is an iteration counter used for the raised 
deck length iteration. The variable moKHI is explained 
1 „ the SKSER subroutine. The t«o values of raised deoh 

length start at zero. 

13 CALL SHEER(KlB.FO.F10,F20.DO.D10.D20,I.EN.H.mCKKT) 

The sheer lino of the vessel is defined by the 
freeboard at stations 0, 10. and 20. These values are 
returned as FO. FiO. and F20 respectively. The associated 

hull depths are also returned. 

AREA«LEN/6 . ♦ (F0+F20+4*F1 0 ) 

FAVG*AREA/i2N 

DAVG«FAVG4H 

The average depth is calculated using the projected 
lateral area. A bettor estimation could be made by using 
the local beam to weight the freeboard values but this 
was not done. because the value of F7 which corrects for 
flare is based on the above average depth. 

F7»-. 01 8828»( DAY G/H)**2+.18098*DAVG/H+. 714599 

This is the flare correction factor. The formula 

a. 1 +hP Navv’s work in developing their 

used was taken from the wavy » 

destroyer model. 

HVAW=LEN*B*CVTP»FAVG*F7 

HVBW*LEN*B*H<»CP*CX 

THV»HVAV/+HVBW 

The total hull volume to the main deck is calculated 
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by those statements* 
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AK-B*H*CX+B*F10 

X«LilB/LEN 

CPM3- , 525*X*X+ . 055*X+ . 99675 
CPM=CPM+3.2«X*X*(CP-.6) 
VMBsLMB*AM*CPM 
AHVOM=THV-VMB 



The volume of the machinery box and the remaining 

available hull volume are calculated in these statei^nts* 

The formula for CPM ha«* been taken from the Navy's work 

and applies to machinery boxes centered about amidships* 

Q*( r , 90371 1 5+ . 21 39?27*D1 0-1 , 38263E-2*D10**2+ 
4.008058B-4i>D10**3-5.489^8lE-6*D10»*4+ 

2. 8921 53E-8^D1 0**5 )<-CP/. 58*CX/* 8l4* 
(1.11*U!B/LEN+,667) 

ATV0M=Q*AHV0M 

VFUEL=WTPUEL*43. 

VL0=WTL0*39. 

VFW=.186*NAC*36. 

VACFUL=V/ACFUL«44, 1 

VREQ* ( VFUEL+YLO+VFW+YACFUL ) *1 . 2 

RTAAaO. 

IF ( VREQ . GT . ATVOM ) RTAA= ( VREQ-ATVOM )/8 . 5 
AAV0K=/iHV'0M-ATV0M 
TANKVL«ATV0?ifRTAA*8 . 5 

The available volume which can be used only for 

tankage eind the required tankage volume are calculated 

and compared in these statements* The formula for Q 

was again taken from the destroyer model work. The 

required tankage arrangements area is calculated if the 

required tankage volume is in excess of that available* 

P* ( 6 * 51 041E-7*D1 0**3-7 . 851484B-5*D1 0**2+3 * 4383753-3* 
D10+* 0530399 )*(.191*CP+*8393)*(.191*CX+* 8445)* 
(-*105*U1B/LBK+1.0315) . . . 

AVAA« AAV 0M*P-RTAA 



The available arrangements volume is converted to 
an available arrangements area again using a Navy developed 
formula* 
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DHV=.00223*LEN**3 

DHAA-DHV/8.5 

CN« ( DAVG+LRD/LEN*8 . 5 ) *LEN*B 
RDHA«( ARSADH+2i| . + . 00225<^N+2, 37*LEN )*1 , 09 
IF(RDHA.GT.DHAA) GOTO 1000 
TAAA=DHAA+AVAA 



The area in the largest acceptable deckhouse is 
compared to the area required to be in the deckhouse.' 

If the required area is smaller, the program continues. 

If not, an error messa^'^ is printed. The areas required 
to be in the deckhouse include the input areas plus the 
commanding officers cabin and the pilothouse. An allowance 
is also made for passageways. 

0DAa2.12*LEN-232. 



PHDA«100,+.00225*CN 

The statements above and those included between 

them are used to calculate the required deck areas. 

TRAA=1 , 09*( AREADH+AREAHL+0DA+J0.IPDA+CSDA+0SRDA+0HDA+ 

C PODA+C POHDA+C BDA+C HDA+C ODA-S-CS TDA+ 03 DA+ WDA+ 
RLDA+CLDA+UTDAtS GDA+ ACDA+XICDA+AMDA+PHDA ) 
IP(TRAA.GT.TAAA) GOTO 10 

The total required arrangements area is computed 

and compared with the area available in the hull to the 

main deck plus the largest acceptable deckhouse. A 

decision is then made to either make the deckhouse smaller 

or to add a raised deck. 

DHSV=. 00151*LEN**3 
DHSAA=DHSV/8. 5 
DHAMIN=AMAX1 ( DHSAA.RDHA) 

TAAA=DHAMI N+ AVAA 
IF(TRAA.GT.TAAA) GOTO 11 
12 DHVaDKSV 
DHAAoDHV/8.5 
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ENCVOLsDHV+THV 

RETURN 

11 DHARsTRAA-AVAA 
DHSV=DHAR*8.5 
GOTO 12 

These statements compute the deckhouse size required 
if the deckhouse is to be made smaller. First, the 
larger value of the minimum deckhouse size or the required 
deckhouse area is chosen and added to the hull voulume. 

This new available area is then compared to the total 
required area. If the available area is still the greater, 
the difference is excess area that cannot be reduced 
since the minimum size hull and deckhouse are already 
being used. If more than the minimum deckhouse size is 
required a calculation is made to balance the total 
required area and the total available area. 

10 RHAA=TRAA-DHAA 
RDAA=RHAA-AVAA 
X1=RDAA/(B*LEH) 

IF(X1.GT. .7) ^OTO 1 
IF(X1.GT..2) GOTO 2 
LRD=-3 . 22*X1 *X1+2 . 031 *X1 - . 0184 
GOTO 3 

1 LRD=.463*X1*X1+.446*X1+.22 
GOTO 3 

2 LRD»X1+.059 

3 LRD=LRD*LEN 
DBCKHT^Se 5 

IF(ABS(LRD-LRD1).LT..05*LEN) GOTO 15 

IF(KNDEX.GT.20) GOTO 1000 

KNDSX»KND2X+1 

LRD1=LRD 

GOTO 13 

If more area is required than is available in the 
largest deckhouse plus the hull to the main deck, a raised 
deck is added. The assumption is made that a ship with 
the maximum deckhouse size will result in the least 







( 



i 



I 

f- 



132 



displacement. This is probably valid so long as an 
aluminum deckhouse is used. 

The length of raised deck required is computed in 
the above statements. The new raised deck length is then 
compared to the former estimate, LRDl, If the two do 
not agree within five percent of the length of the ship 
a new iteration is made A default counter has been 
included to keep the computer from being caught in a 
nonconverging loop. 

15 C1«.4*LBN 
C2=.6*LEN 

IF(LRD.LT.C1.0R.LRD.GT.C2) GOTO 4 
AVAA=AVAA+ ( C2-LRD) 

RDAA*RDAA+ ( C 2 -LRD ) 

LRD=C2 

These statements check to see if the raised dock 
length falls within 0,4 to 0,6 times the ship length. 

If so the calculated values are corrected so that the 
raised deck length is equal to 0,6 times the ship length, 

4 IF(LRD.GT.LEN) GOTO 1000 
THV=THV+RDAA«8 , 5 
ENCV0L»DHV+THV 
RETURN 

A check is first made to insxure that a raised deck 

length longer than the length of the ship is not required, 

A corrected total volume is then calculated and the program 

returns, 

1000 LEN=-1000. 

RETURN 

END 

Statement 1000 is a default option warning subroutine 
XECUTE to print 2 tn error message. 
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3.10.^ Output List for Subroutine VOLUME 



ACDA. . . . 




DHV 




AMDA. . . . 




ENCYOL 




AVAA.... 




LRD 




CBDA. . . . 




ODA. 




CHDA.... 




OHDA 




CLDA. . . . 




OSDA 




CN 




OSRDA 




CODA.... 




PHDA 




CPODA... 




RLDA 




CPOHDA. . 




SGDA 




CSDA. . . . 




TANKVL 




CSTDA. . . 




UTDA 




DO 




VHB 




DIO 


LC(II ) 


WDA 


LC(PP) 


D20 




XICDA 




DAYG. . . . 




XMPDA 




DHAA. .. . 









3*10.5 Nonenclaturo List 

The definition of all variables is the same as 
given in the MAIN program nomenclature except for the 
following! 



AAVOM 

ACDA 

AHVOM 

AM 

AMDA 

AREA 



available arrangements volume in hull, cu ft 
A/C & ventilation deck area, sq ft 
available hull volume , cu ft 
cross sectional area amidships, sq ft 
auz. machinery spaces deck area, sq ft 
projected lateral area above water, sq ft 
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ATVOM 

AVAA 

Cl 

C2 

CBDA 

CHDA 

CLDA 

CODA 

CPM 

CPODA 

CPOHDA 

CSDA 

CSTDA 

D 

DO 

DIO 

D20 

DECKHT 

DHAA 

DHAMIN 

DHAR 

DHSAA 

DHSY 

1MB 

FO 

F? 

FIO 



available tankage voluse, cu ft 

available arrangements area in hull, oq ft 

0.4 X length, ft 

0,6 X length, ft 

crews berthing deck area, sq ft 

crews heads deck area, sq ft 

cl^in locker deck area, sq ft 

C.O. S.R., Cabin & pantry deck area, sq ft 

machinery box prismatic coefficient 

CPO S.R. deck area, sq ft 

CPO heads deck area, sq ft 

crew special deck area, sq ft 

commisoary stores deck area, sq ft 

same as ENDDAY in MAIN 

depth at F.P. , ft 

depth at amidships, ft 

depth at A.P. , ft 

height of raised dock, ft 

deckhouse arrangements area, sq ft 

minimum deckhouse area, sq ft 

deckhouse area required, sq ft 

smallest deckhouse arrangements area, sq ft 

smallest deckhouse volume, cu ft 

same as DHM in MAIN 

freeboard at F.P., ft 

flare factor 

freeboard at amidships , ft 
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F20 


freeboard at A.P. ,ft 


FAVG 


average freeboard, ft 


HVAW 


hull volume above waterline, cu ft 


HVBW 


hull volume below waterline, cu ft 


KNDEX 


counter 


LMB 


same as XLM in MAIN 


LRDl 


used to store LRD value 


NAG 


number of accommodations. 


ODA 


office deck area, sq ft 


OHDA 


officer heads deck area, sq ft 


OSDA 


other stores deck area, sq ft 


OSRDA 


officer staterooms deck area, sq ft 


P 


coefficient 


PHDA 


pilothouse, chartroom A CIC deck area, oq ft 


Q 


coefficient 


RDAA 


raised dock arrangements area, sq ft 


RDHA 


area required to be in deckhouse , sq ft 


RE 


same as RGBND in MAIN 


RHAA 


required hull arrangements area, sq ft 


RLDA 


repair locker deck area, sq ft 


RTAA 


required tankage arrangements area, sq ft 


SFCH 


same as SFCHHP in MAIN 


SFCM 


same as SFCMHP in MIN 


SGDA 


steering gear & windlass deck area, sq ft 


TAAA 


total available arrangements area, sq ft 


TANKVL 


tankage volume , cu ft 
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THV 


total hull voluiae , cu ft 




TRAA 


total roquirod arrangemonts 


aireat 


UTDA 


uptake deck area in hull, sq ft 


VACFUL 


volume of aircraft fuel, cu 


ft 


VE 


same as VEND in MAIN 




VFUEL 


volume of ships fuel, cu ft 




VFW 


volume potable water, cu ft 




YLO 


volume of lub oil, cu ft 




VMB 


volume of machinery box, cu 


ft 


YREQ 


required tankage volume , cu 


ft 


WDA 


workshops deck area, sq ft 




X 


LMB/LEN 




XI 


RDAA/(B«LEN) 




XICDA 


I.C, & gyro room deck area. 


sq ft 


XMFDA 


messing facilities deck area, sq : 
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3.11 Subroutine SHEER 
3 , 11,1 Introduction 

Subroutine SHEER selects an appropriate sheer line 
by specifying the freeboard at the forward perpendicular, 
amidships, and at the after perpondic\aar. This sheer 
line must satisfy a nxxmber of criteria. 

First, there must be sufficient depth amidships 
to accommodate the engine room. The depth amidships 
must also be sufficient to insure adequate structural 
strength, A minimum veiluo of l/i 6 is used. Because 
Coast Guard vessels normally have engine rooms aft of 
amidships, a slightly greater depth amidships is used 
rather than the machinery box depth. This accounts for 
the drop in the sheer line aft of amidships. 

The next criterion which must be satisfied is adequate 
freeboard at the forward and after perpendiculars. 

Navy derived formulas were used. These formulas give 
the minimum acceptable freeboard based on a dock wetness 
criterion. 

Once independent determinations of freeboard at the 
forward and after perpendiculars and at aiaidships are 
made , a check is required to insure that a reasonable 
sheer line can be drawn between the three. Data was 
obtained from past designs and from the standard merchant 
ship sheer curves. This data was nondiraonsionalized by 
subtracting the freeboard amidships from the other free- 
board values amd then dividing by the length of the ship. 
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The resulting curves are plotted in Figure 28, 

The sheer lino is assumed to be practical if the 
forward sheer fraction lies between 0,01 and 0,03 and the 
after sheer fraction lies between 0,001 and 0,0075* 

If this is not the case, one or more of the freeboard 
values is adjusted until the criterion is satisfied. 

If a raised deck is added to the hull, a new sheer 
line is calculated. It is assumed that the required 
freeboard of the main deck at the forward perpendicular 
can bo one deck height or 8,5 feet below that required 
by the deck wetness criterion. This will result in a 
flattened sheer lino as is common for raised deck vessels, 
A flow chart for this routine is shown in Figure 29* 

3,11,2 Input List for Subroutine SHEER 

DMB SCDA H SCDA 

LEN SCDA DECKHT SCDA 

3. 11 *3 Statement Descriptions 

D10o,0025*LEN+K.lB 
IF(D10.LT.LEH/18) D10 bLEN/16, 

The minimum depth amidships is selected by these 
statements, 

F0=l, Oil 827* (1 00, *H/LEN)-,000636215*LEN+2, 780649 
F0=LEK«F0/100, 

FOaFO-DSCKHT 

DOaPO+H 

The minimum depth at the forward perpendicular is 



selected here. 
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F20-.01*LEN*(2.125+.00125«^LEN) 

D20oF20+H 

Finally the nininmin depth at the after i>erpendicular 

is selected. 

7 SOaDO-DlO 
R0=30/LEN 

IF(RO.LT..Ol) GOTO 1 
IF(R0.GT..03) GOTO 2 

The sheer fraction forward is calculated and checked 

against the acceptable limits, 

6 S20=D20-D10 
R20=S20/LSN 

IF(R20.LT. .001) GOTO 3 
IP(R20,GT. .0075) goto k 
GOTO 5 

The sheer fraction at the after perpendicular is 

calculated and checked against acceptable limits. 

1 S0=.01*LEIi 
DO-SG+DIO 
GOTO 6 

If the freeboard at the bow is too low the depth 
at the bow is increased. 



2 S0s.03*LEN 
D10=D0-S0 
GOTO 6 

If the freeboard at the bow is excessive , the depth 
amidships is increased, 

3 S20a.001*LEN 
D20=D10+S20 
GOTO 5 

If the freeboard aft is too low the depth aft is 
increased. 

4 S20=.0075*LEN 
D10=D20-S20 
GOTO 7 
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If the froehoard aft is excessive the depth amidships 
is increased and the sheer fraction at the bow is rechecked 



5 FOsDO-H 
F20=D20-H 
FlOsDlO-H 
RETURN 
END 

Values of both freeboard and depth are returned 
for use in subroutine v^olume. 

3,11.4 Output List for Subroutine SHEER 



PO 


SODA 


DO 


...SODA 


FIO 




DIO 





F20 SODA D20 SODA 



3.11.5 


Nomenclature List 


DO 


depth at forward perpendicular, ft 


DIO 


depth at amidships, ft 


D20 


depth at after perpendicular, ft 


DECKHT 


height of raised deck, ft 


DC4B 


machinery box depth, ft 


PO 


freeboard at bow, ft 


FIO 


freeboard amidships, ft 


F20 


freeboard at stem, ft 


H 


dreift, ft 


LEN 


length between i>erpendiculara , ft 


RO 


(D0-D10)/LEN 


R20 


(D20-D10)/I£N 


SO 


DO-DIO 



S20 



D20-D10 
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3.12 Subroutine WEIGHT 
3,12,1 Introduction 

Subroutine V0LUI*1E balanced the required and available 
volumes in the hull. This subroutine estimates weights 
so that a balance can be made between full load weight 
and displacement. 

Several of the weights required have already been 
calculated or were given as input. Calculated weights 
include the fuel weight and lub oil weight. All com- 
munications and control, group 4, weights and all armament, 
group 7, weights are given as input. The weight groups 
noted refer to the NAVSHIPS Hull Group Weight Classification 
of 1965. The categories in this classification scheme 
are given in Appendix E. This is the classification 
system now used by the Coast Guard. 

Also given as input are the weights of cargo, air- 
craft, ammunition and aircraft fuel. All other weights 
which are a part of the full load weight are estimated 
in subroutine WEIGHT. A flow chart for this routine is 
.shown in Figure 30* 

There are few patrol cutters of recent vintage 
which can be used for collecting weight data. Three 
have been used throughout this model. These are the 
378’ WHEC, the 210’ WMBC, and the 150* WPB. Wherever 
possible additional data has been shown. Light ship 
weight data was taken from the weight breakdown sheets 
prepared for each design. Load weights were taken from 
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either load weight calculations or froa inclining experiaont 
data. 

The first hull weight group is the hull structure, 
group 1» In order to estimate this freight group more 
accurately, three of the sub-groups are estimated sep- 
arately. The remainder of the sub-groups are lumped 
together as a single group. In the program this group 
is referred to as WTIH. The three groups subtracted 
are group 111, superstructure! group 112, foundations 
for propulsion plant machinery! and group 113 » foundations 
for auxiliaries and other equipment. These are identified 
in the prograira as WTlll, WT112, and WT113, respectively. 

Figure 31 is a plot of WTIM versus cubic number. 

The correlation hero is very good. Also included on the 
figure is the estimating lino used in tho Navy's destroyer 
model, DDO 7 . 

Superstructure weight, WTlll, is shown in Figure 32. 
plotted against deckhouse volume. The line labeled no 
uptakes is the one used in the program. The line used 
in DDO 7 is also shown. This line assumes that the deck- 
house is aluminum. 

The weight of propulsion plant foundations, WT112, 
is shown in Figure 33 plotted against the weight of 
propulsion machinery minus shafting and propellers. Again, 
a line used in DD 07 is shown, 

WTII 3 , other foundations, is shown in Figure 3^ 
versus the total weight of groups 3» 5» and 7, The 
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correlation here appears reasonable; howeveri there is 
a difference betreen the data and the DDO7 estimate. 

Propulsion machinery is classified under group 2. 

In this program, an estimate of the weight of sub-group 
203 (shafting, bearings, and propellers) was made separately. 
The remaining sub-groups were estimated as one group, WT2M. 

A plot of WT2M versus shaft horsepower for diesel plants 
is shovm in Figure 35 • Here more data is available from 
non patrol types. The DD07 estimate for diesel plants 
is poor. For CODOG plants it is assumed that the diesel 
plant will be 7OOO horsepower or under and that only 
the endurance power will cause a difference in weight 
from that of a 378’ WHBC. A lower endurance power will 
reduce the weight somewhat. 

Shafting and propeller weight, WT203, is calculated 
as weight per foot of shaft length versus sheift horsepower. 
Twin shafts are assumed vrithatotal shaft length O.56 
times the length of the ship. A plot of this relationship 
is shown in Figure 36. 

Group 3 includes the electrical plant weights. The 
Coast Gusird intends to install electric plants of relative- 
ly leirger capacity in future ships so a simple linear 
estimation using current designs is unrealistic. Instead, 
the weights of a number of diesel generating sets of 
various capacities have been plotted as shown in Figure 37* 
The varying RPM of those plants resulted in three weight 
ranges. 
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The weight of the generating sets is only a part 
of the group 3 weight. The weight of switchboards, WT301, 
is shown in Figure 38 as is the remaining group 3 weight, 
WT3R. The correlation for those weight groups is quite 
bad. Fortunately, these weights make up a small percentage 
of the total weight of the ship. 

Since group 4 weignts are given as input, group 5 
weights, auxiliary systems, is estimated next. The 
destroyer model uses the estimating relation shown in 
Figure 39. The three Coast Guard vessels agree very 
well with this relationship . 

Group 6, outfit and furnishings is the last light 
ship group estimated since group 7 weights are given as 
input. Group 6 is broken into four sub-groups for estimating 
purposes. The fii*st, WT61, includes groups 6OO, 602, 

603, and 605. These are hull fittings, rigging and canvas, 
ladders and gratings, and painting, respectively. The 
second group, WT611, includes groups 604, nonstruotural 
bulkheads and doors 1 606, deck coveringi 607, hull insulation 
and 608, storerooms, stowages and lockers. Group three, 
WT6111, includes only group 601, boats, boat stowage 
and handling. The final group, WT61V, includes all the 
remaining sub-groups of weight group 6. 

WT6i plotted against enclosed voltine is shown in 
Figure 4o. WT611 and WT6IV are plotted against number of 

accommodations as shown in Figure 41. WT6111 is estimated 

in three parts. The first is a boat weight, the second a 
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liferaft weight and the third a lifejacket weight. 

All of the lightship weights are sumaed and then 
multiplied by a design and builders margin. The design 
and builders margin is an input. 

Crew weight is taken to be 0.0737 tons per man or 
165 pounds per man. Personeil effects are calculated 
as 0,105 tons per man for officers, .0737 for chiefs 
and .039 for enlisted men. 

The weight of personnel stores includes provisions, 
general stores, and potable water. The weight of provisions 
is based on Table 3* The data in this table was reduced 
to a formula for which only the endurance days of dry 
provisions need bo given. This formula isi 
(.0222+. 00202 x D) x NAG 

where D is the enduraince days for dry provisions 
NAC is the number of accommodations 

Table 3 

days endurance pounds/person 





dry prov. 


frozen 


chilled 


total 


30 


97.0 


36.8 


51.8 


185.6 


45 


145.5 


55.2 


52.8 


253.5 


60 


194.0 


73.6 


53.8 


321.4 


90 


291.0 


110.4 


55.8 


457.2 


The 


weight of general 


stores is 


taken to be 





.00135 X D X NAC 
and that of potable water to be 



.186 X NAC or 50 gallons/man 
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The full load vreight can then be calculated by adding 
the load weights to the light ship weight. 

3.12,2 Input List for Subroutine WEIGHT 



CN 




SHPB 




D 




SHKl 




DBMAR... 




WACFUL 




DHV 




WCARGO 




ELKW.... 




WTAC 




ENCVOL.. 




WTAKMO 




JOPT. . . . 




WTFUSL 




LEN 




WTG2T 




MTyPE... 




WTG4, 




NCPO.... 




WTG7 




NENL.... 




WTLO 




HOFF. . . , 










♦only if JOPT «. 2 







3.12,3 Statement Descriptions 

A detailed description of each statement will not 

be given for this routine wince the flow of the program 

follows the description given in the introduction and 

the relationships are those plotted. 

One area that will be discussed is WT6111. Here 

the nximber of 15 man liferafts required is first calculated, 

Q*NAC/15.+1. 

RAFTSaAINT(Q) 

where AINT is a function which drops the decimal fraction 
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from its argument. Boat weight is assumed to be 1.4 tons 
up to 25 men, 5 tons between 25 and $0 men and 10.1 tons 
above $0 men. This weight includes launching equipment* 
Lifejacket weight is taken to be .002 tons per man. 

3.12.4 Output List for Subroutine WEIGHT 



WFULLD 




WT 611 




WLSHIP 




WT 6 III 




WTlll 




WT 6 IV 




WT112 




WTCREW. . , . . . 




WTII 3 




WTGl 




WTIM 




WTG2T 




WT203 




WTG3 




WT2M 




WTG5 




WT 3 OO 




WTG6 




WT 3 OI 




WTPE 




WT3R 




WTPS 




WT 6 I 









♦ only if JOPT = 1 



3 . 12.5 Nomenclature List 

The definition of all variables is the same as 
given in the MAIN program nomenclature except for the 
foil owing I 

BOATWT boat weight, tons 

D same as ENDDAY in MAIN 

NAC number of accommodations 

dummy argument 
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RAFTS 



RE 

SFCH 

SFCM 

VE 

VOL 

WTlll 

WT112 

WTII3 

WTIM 

WT203 

WT2M 

WT300 

WT301 

WT3R 

WT6i 

WT 6 II 

WT 6 III 

WT6IV 

WTFDS 



number of liferafts 

same as RGEND in MAIN 

same as SPCHHP in MAIN 

same as SFCMHP in MAIN 

same as VEND in MAIN 

enclosed volume x 10”^ 

weight of group 111, tons 

weight of group 112, tons 

weight of group 113 » tons 

weight of remainder of group 1, tons 

weight of group 203 , tons 

weight of remainder of group 2, tons 

weight of group 300, tons 

weight of group 30 I, tons 

weight of remainder of group 3, tons 

see text 

see text 

800 text 

see text 

sum of weight groups 3 , 4, 5 , and 7 
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3.13 Subroutine VCO 
3,13.1 Introduction 

The vertical center of gravity of the ship is estimated 
in this routine. The designer has some latitude with 
regard to the location of the center of gravity; there- 
fore, the value calculated by this routine should be viewed 
as a feasible rather tnan as the best value. 

The effect that the veg has on the design can be 
investigated by using different values of free surface 
correction. The free surface correction is added to the 
KG. Therefore, either a positive or a negative value 
can be used to represent a shift in the KG location. 

In order to determine feasible locations for the 
various centers, data for the 378' WHEC, the 210* V/MEC, 
and the 15O' WPB was used. This data is shown in Tables 
4 and 5» On these tables the vegs of the various weight 
groups are shown. These centers are related to various 
hull depths by a factor K. The hull depth used is given 
in the right hand column of these tables. An average 
coefficient, K, was chosen and is given in the column 
labeled Coef, K. 

The variables to which the vegs are referred were 
chosen arbitrarily by observing which gave the most 
consistent set of coefficients among the three ships. 

Some of the weight groups shown in the tables are 
really input items. These weight groups are included 
to assist the designer in choosing realistic values to 
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Weight 




373' 

WHEC 


210* 

WIvEC 


150' 

WPB 


Coef . 
K 


times 


WTli'i 


vcg 

K 


15.73 

.553 


14.72 

.547 


9.55 

.584 


.55 


DAVG 


WTlll 


vcg 

K 


37.1 

1.4 


35.3 

1.26 


23.6 

1.47 


1.4 


DIO 


WT112 


vcg 

K 


5.73 

.216 


4.4 

.157 


2.5 

.156 


.2 


DIO 


WT113 


vcg 

K 


21.05 

.794 


16.6 

.593 


15.1 

.944 


.87 


Center groups 
3+4+ 5+ 7 


WT2I.: 


vcg 

K 


15.14 

.571 


10.40 

.533 


8.50 

.531 


.54 


DIOMN 


WT203 


vcg 

K 


4.02 

.152 


6.0 

.214 


2.6 

.163 


.16 


DIO 


WT300 


vcg 

K 


20.35 

.768 


14.7 

.525 


19.0 

1.19 


.7 


DIO 


WT301 


vcg 

K 


21.3 

.804 


11.8 

.421 


21.0 

1.31 


.8 


DIO 


WT3R 


vcg 

K 


27.5 

1.038 


19.41 

.693 


17.79 

1.11 


1.0 


DIO 


WTG5 


vcg 

K 


19.43 

.733 


1 6 • 63 
.596 


13.5 

.844 


.71 


DIO 


WToI 


vcg 

K 


20.66 

.78 


20.43 

.73 


15.09 

.943 


.75 


DIO 


WToII 


vcg 

K 


26.48 

1.0 


21.93 

.783 


19.23 

1.2 


1.0 


DIO 


wt6iii 


vcg 

K 


41.75 

1.575 


33.9 

1.738 


25.9 

1.62 


1.6 


DIOIVIN . 


wt6iv 


vcg 

K 


24.88 

.939 


20.97 

.749 


13.46 

.841 


.84 


DIO 


■WTG4 


vcg 

K 


31.56 

1.191 


35.0 

1.25, 


22.6 

1.41 


1.3 


DIO 


V/TG? 


vcg 

K 


35.26 

1.33 


24.45 

.873 


21.3 

1.33 


1.3 


DIO 



Table 4 
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Weight 




378’ 

WHSC 


210’ 

WIvIEC 


150* 

WPB 


Coef . 
K 


times 


Crew & 

personal 

Effects 


vcg 

K 


26.26 

.964 


20.43 

1.048 


15.0 

.938 


.97 


DIOMN 


Provisions vcg 
K 


15.0 

.55 


15.0 

.536 


11.0 

.688 


.55 


DIO 


General 

Stores 


vcg 

K 


24.05 

.883 


16. 16 
.829 


11.0 

.688 


.85 


DIOMN 


Ammo 


vcg 

K 


15.99 

.587 


8.4 

.431 


19.0 

1.19 


.5 


DIOMN 


Potable 

Water 


vcg 

K 


7.06 

.26 


8.5 

.30 


4.2 

.263 


.26 


DIO 


L.C. 


vcg 

K 


3.58 

.131 


2.5 

.128 


12.0 

.75 


.13 


DIOMN 


Fuel 


vcg 

K 


7.67 

.281 


3.17 

.292 


6.6 

.412 


.3 


DIO 


Aircraft 


vcg 

K 


39.01 

1.43 






1.4 


DIO 


A/C Fuel 


vcg 


9.56 

.351 






.351 


DIO 


DAVG 




28.54 


26.91 


16.36 






DIO 




27.25 


28.0 


16.0 






D1 0Iv:N 




27.25 


19.5 


16.0 







Table 5 
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Figure 42 b 
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input 
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l^ZsL and 42b. 

3.13.2 

ACCG 
ACWT 
AMXG 
AMOWT 
CACFUL 
CARGOC 
CARGOW 

D 

DIOMN. 

DAVG. . 

DBMAR. 

GR4CG. 

GR4WT. 

GR7CG. 

GR7WT. 

LEN... 

LRD. . . 

NOARM. 

NOELT. 

WACFUL LC(LL) 

WCARGO LC(LL) 

WFULLD LC(WW) 



mi 3 




WTIM 




WT203 




WT2M 




WT300 




WT 3 OI 




WT3R 




WT 6 I 




WT 61 I 




WT 6 III 




wt 6 iv 




WTAMMO 




WTCREW 




WTPUEL 




WTGl 




WTG2 




WTG3 




WTG4 




WTG5 




WTG 6 




WTG7 




WTLO 





A flow chart for this routine is shown in Figures 
42b. 

Input List for Subroutine VCG 
LC(EE) 

LC(EE) 

LC(EE) 

LC(EE) 

LC(HH) 

LC(EE) 

LC(EE) 

LC(AA) 

LC(II) 

LC(BB) 

LC(WV/) 

LC(EE) 

LC(EE) 

LC(EE) 

LC(EE) 

LC(BB) 

LC(JJ) 

LC(FF) 

LC(FF) 



WLSHIP 




WTPE 
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WTlll 




WTPS 




WT112 








3.13.3 Statement Descriptions 







Detailed statement descriptions will not be given for 
this routine. In general the statements consist of vcg 
estimates using the relationships mentioned in the intro- 
duction or of calculations to combine weight groups by 
summing the vertical moments of individual groups and 
dividing by the total weight. 

Note that DIO is different from DIOMN only if the 
length of raised deck is greater than one half the ship 
length. It must also be remembered that input moment 
arms are given as fractions of the main deck depth at 
amidships and not as a distance in feet. 



3.13.4 Output List for Subroutine VCG 



CCARGO LC(OO) CGG3 

CFULD SODA CGG4 

CFULLD LC(OO) CGG5 

CGAC LC(HH) CGG6 

CGAI^MO LC(HH) CGG? 

CGCREW LC(HH) CGLO 

CGFUEL LC(HH) CGPE 

CGGl LC(OO) CGPS 

LC(OO) 



LC(OO) 

LC(OO) 

LC(OO) 

LC(OO) 

LC(OO) 

LC(HH) 

LC(HH) 

LC(HH) 

LC(OO) 



CGG2 



CLSHIP 
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3.13.5 

ACMOM 

AMOMOM 

CACFL 

CFULD 

CFULLD 

CGlll 

CG112 

CGII 3 

CGIM 

CG203 

CG2M 

CG300 

CG 3 OI 

CG3R 

CG6I 

cg6ii 

CG 6 III 

CG 6 IV 

CGGl 

CGG2 

CGG3 

CGG4 

CGG5 

CGG6 

CGG7 



total vertical aircraft moment, tons 
total vertical ammunition moment, tons 
vortical center aircraft fuel, ft 
same as CFULLD 



vertical 


center 


of 


ship at full load. 


vertical 


center 


of 


group 


111 


. ft 


vertical 


center 


of 


group 


112 


. ft 


vertical 


center 


of 


group 


113 


. ft 


vertical 


center 


of 


group 


IM, 


ft 


vortical 


center 


of 


group 


203 


. ft 


vertical 


center 


of 


group 


2M, 


ft 


vertical 


center 


of 


group 


300 


. ft 


vertical 


center 


of 


group 


301 


. ft 


vertical 


center 


of 


group 


3R. 


ft 


vertical 


center 


of 


group 


61 . 


ft 


vertical 


center 


of 


group 


611 


. ft 


vertical 


center 


of 


group 


6111 , ft 


vertical 


center 


of 


group 


61 V 


r ft 


vertical 


center 


of 


group 


1. 


ft 


vertical 


center 


of 


group 


2, 


ft 


vertical 


center 


of 


group 


3. 


ft 


vertical 


center 


of 


group 


4, 


ft 


vertical 


center 


of 


group 


5. 


ft 


vertical 


center 


of 


group 


6 . 


ft 


vortical 


center 


of 


group 


7. 


ft 


vertical 


center 


of 


general stores, ft 



CGGS 



1 

\ 



I 
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CGOMOH 


verticeQ moment of cargo, tons 


COPY 


vertical center of provisions, ft 


CGPW 


vertical center of potable water, ft 


CLSHIP 


vertical center of light ship, ft 


CPO 


same as NCPO in MAIN 


D 


same as ENDDAY in MAIN 


DO 


depth at F.P., ft 


DIO 


depth amidships to main or raised deck, ft 


DIOMN 


depth amidships to main deck, ft 


D20 


depth at A.P. , ft 


ENL 


same as NENL in MAIN 


G4K0M 


total group 4 moment, tons 


G7MOM 


total group 7 moment, tons 


HLEN 


half the length of the ship, ft 


I 


index 


OFF 


same as NOFF in MAIN 


WTIM 


weight group IM, tons 


WTlll 


weight group 111 , tons 


WT112 


weight group 112 , tons 


.WTII 3 


weight group 113 » tons 


WT203 


weight group 203 , tons 


WT2M 


weight group 2M, tons 


WT300 


weight group 300 , tons 


WT 3 OI 


weight group 301 » tons 


WT3R 


weight group 3 R» tons 


WT 6 I 


weight group 6l , tons 


WT 61 I 


weight group 611 , tons 



I 

{ 



i 



WT6III 

WT6IV 

XLDMOM 
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weight group 6lII , tons 
weight group 6lV, tons 
total load moment, ton-ft 
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3.14 Subroutine COST 
3.14,1 Introduction 

This subroutine is based exclusively on the pre- 
liminary cost estimating procedure currently used by the 
Coast Guard, To use the procedure, the weight of each 
of the seven light ship weight groups must be known. 

Linear graphs for man hours and material cost are then 
entered with these values to determine the cost and labor 
required for each weight group. 

Various relationships are also given for determining 
design costs, construction services cost and miscellaneous 
costs. Miscellaneous costs include initial outfit equip- 
ment such as lifejackets as well as retrofit costs, spare 
parts and administrative or resident inspector’s office 
costs. 

The curves used in the procedure were converted 
to formulas and the procedure was programined as subroutine 
COST. A flow chart for this routine is given in Figure 

43. 

All of the material cost curves are based on 1959 
prices. Indices are used to correct the prices to the 
year the ship will be built. These indices are based 
on Bureau of Labor Statistics data for the particular 
weight group. 

The labor rate for all of the weight groups is assumed 
to be the same. The labor rate is inputed as an average 
dollars per hour for the year the ship will be built. 



1 

1 

i 

I 

i 

I 



COST SUBROUTINE 
FLOW CHART 
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Material cost indices and labor rate data are shown 
in Figures 44 and 45, There is a large fluctuation in 
costs between shipyards so these numbers can only be 
viewed as average values. 

In addition to the corrections for inflation, eighty 
percent is added to labor cost for overhead and ten 
percent is added to all costs for profit. 

Initial outfit costs are assumed to be two percent 
of the total cost, Spaure parts cost is the sum of the 
cost of spares for propulsion machinery, electrical 
machinery, electronics and outfit. The cost of spares 
as a percentage of the weight group material cost is 
given be low I 

Group 2 9 % 

Group 3 — 8 % 

Group 4 25^ 

Group 5 - — 10^ 

Retrofit cost is assumed to be 4 percent of the 
total cost. Administrative socts are assumed to be 3,5 
percent of the total. 



3.14,2 Input List for Subroutine COST 



LC(MM) 


GlIND 




LC(M>!) 


G2IND 




LC(MM) 


G3IND 




LC(WW) 


G5IND 




LC(NN) 


G6IND 





DOLHR 



i 

I 



INDEX 
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YEAR 
Figure 44 



5.00 
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dH /$ 31VU uoavi 



YEAR 
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JOPT 




WTG3 




JJEU 


.a .SODA 


WTG4. . . . . 


...... .LCfWVf) 


MTYPE 




WTG5 




WTGl 




WTG6 




WTG2 




WTG7 





♦only if JOPT « 2 



3. 1^*3 Statement Descriptions 
W1=WTG1*DBMAR 

W7=WTG?*DBMAR 

The design and huilders margin is evenly divided 
among the weight groups for costing purposes. 
G1MH=10.**(AL0G10(W1)*. 75217+2. 97276) 

• 

G7Wi=l 0 .**( ALOGl 0 ( W7 )*. 7521 7+2 . 97276 ) 

These formulas represent the relationships between 
man-hours and the weight of individual weight groups. 
These relationships are assumed linear when plotted on 
log-log paper. 

G1MC=10.**(AL0G10(W1 )♦. 95^243+2. 79 0^8if) 

G2MC=CSTG2 

IF(J0PT.EQ.2) GOTO 1 

G2KC=10.**(AL0G10(W2)*1. 018885+3.634328) 
IP(MTYPE.EQ.2) G2MC«G2MC+1500000. 

1 G3MC=10.**(AL0G10(W3)*1. 072551+3.667812) 

G4MC=CSTG4 

G5MC«10.**(AL0G10(W5)*. 870906+3. 974192) 

G6MC=1 0 .*♦( ALOGl 0 ( W6 ) *1 . 068265+3 . 603833 ) 

G7yiC*CSTG7 

These formulas represent the material cost curves 
used by the Coast Guard. Some of these costs. CSTG2, 
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CSTG4 and CSTG7, are input. Group 2 cost is an input 
only if the machinery plant is spocifiod. Otherwise, 
the machinery cost is calculated. The fonnula used 
represents the cost of diesel plants. An extra $1.5 
million is added if a CODOG plant is to be installed. 
G1LC=G1MH*D0LHR*1 . 98 

G7LC=G7MH»DOLHR«l . 98 

The labor cost by weight group is calculated here. 

The value, 1.98, represents an overhead factor of 1,8 

times a profit factor of 1.1. 

G1KC=G1MC«G1IND»1 . 1 
IF(J0PT.EQ.2) GOTO 2 
G2KC»G2MC*G2IND 
2 G3MC=G3r‘'C*G3IND*l.l 
G5MC=G5MC *G5IND*1 . 1 
G6MC=G6MC*G6IND*1 . 1 
G2MC=G2MC*1.1 
G^MC*G4MC*1.1 
G7MC=G7?'IC*l.l 

In these formulas the material cost is corrected 

for inflation and a factor of 1.1 is included for profit. 

Note that the costs given as input are not corrected 

for inflation. Therefore, the cost of the equipment 

must be corrected for inflation before being inputed. 

LABCST=G1LC+G2LC+G3LC+G4lC+G5LC+G6lC+G7LC 

MATCST=GlMC+G2MC+G3aC+G4KC+G5MC+G6MC+G7MC 

The totad labor and material costs are calculated 



in these statements 
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L0A=1.08*LEN 

FACTORS (.114*L0A-1, 6 )/l 00. 

DLC =? AC TOR *L ABCS T 
FACTOR=( - . 01 «L0A+7 . 5 )/lOO. 

DMCsFAC T0R*MATCS T 

These statements calculate the design labor, DLC, 

and design materials, DMC, costs for the first ship of 

a class from the beginning of the conceptual design stage. 

The two values of FACTOR are based on Coast Guard curves. 

Since the length overall is required in these curves, its 

value is estimated as eight percent greater than the 

length between perpendiculars, 

FACTOR* (.014*L0A+7. )/100, 

CLC*FACTOR«-LABCST 

FACT0R=(-.015*L0A+6.5)/100. 

CKC»FACTOR«MATCST 

The construction services labor and material costs 

are calculated here. Again the formulas for the two values 

of FACTOR are based on Coast Guard curves. 

TOTCS TsL/iBCS T+MATCST+DLC+DI.1C+CLC+CMC 
0RICST=.02»T0TCST 

SPARES* ( . 09*G2MC+ . 08«G3MC+. 25<^G4MC+. 1«G5MC )/l , 1 
RETRO*. 04»T0TCST 
ADMIN*. 035*T0TCST 

The four items that are included in tho miscellaneous 
cost are calculated here. Since profit is included in 
the individual weight group material costs it is divided 
out when calculating the cost of spares, 

CSTGI-GILC+GIHC 



CSTG7-G7LC+G7MC 

Total weight group costs are calculated here 






I 



i 
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DCST=DLC+DMC 

CCST=CLC+CMC 

MCSTaORICST+SPARES+RETRO+AKdIN 

TOTCST»TOTCST+MCST 

RETURN 

END 

Finally, the total design cost, construction services 
cost, miscellaneous cost and the total first ship cost are 
calculated* 



3*1^.^ Output List for Subroutine COST 



CCST 




CSTG6 




CSTGl 




CSTG7 




CSTG2t 




DCST 




CSTG3 




MOST 




CSTG4 




TOTCST 




CSTG5 









♦ only if JOPT = 1 



3*14*5 Nomenclature List 

The definition of all variables is the same as that 
given in the MAIN program nomenclature except for the 
following I 



ADT«N 

CLC 

CMC 

DLC 

DMC 

FACTOR 

GILC 



administrative costa, dollars 
construction services labor cost, dollars 
construction services material cost, dollars 
design labor cost, dollars 
design materials cost, dollars 
see text 

weight group 1 labor cost, dollars 






4 

4 

I 

( 
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G2LC 

G3LC 

G4LC 

G5LC 

G6LC 

G7LC 

GIMC 

G2MC 

G3MC 

G4MC 

G5MC 

G6MC 

G7MC 

GIMH 

G2MH 

G3MH 

G4MH 

G5MH 

G6MH 

G7MH 

LABCST 

LOA 

MATCST 

ORICST 

RE 

RETRO 

SFCH 



weight 


group 


2 


labor cost, 


dollars 


weight 


group 


3 


labor cost. 


dollem^s 


weight 


group 


4 


labor cost. 


dollars 


weight 


group 


5 


labor cost. 


dollars 


weight 


group 


6 


labor cost. 


dollars 


weight 


group 


7 


labor cost. 


dollars 


weight 


group 


1 


material cost. 


dollars 


weight 


group 


2 


material cost. 


dollars 


weight 


group 


3 


material cost. 


dollars 


weight 


group 


4 


material cost. 


dollars 


weight 


group 


5 


material cost. 


dollars 


weight 


group 


6 


material cost. 


dollars 


weight 


group 


7 


material cost. 


dollars 


weight 


group 


1 


man-hours 






weight 


group 


2 


man-hours 






weight 


group 


3 


man-hours 






weight 


group 


4 


man-hours 






weight 


group 


5 


man-hours 






weight 


group 


6 


man-hours 






weight 


group 


7 


man-houirs 







total construction labor cost, dollars 
length overall, ft 

total construction material cost, dollars 

initial outfit cost, dollars 

same as RGEND in MAIN 

retrofit cost, dollars 

same as SFCHHP in MAIN 



i 

{ 

I 

( 

( 

i 

( 

i 

4 

i 
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SFCM 

SPARES 

YE 

W1 

W2 

W3 

W4 

W5 

W6 

V7 



same as SPCMHP in MAIN 
spare parts cost, dollars 
same as VEND in MAIN 
vreight group 1 with margin, 
weight group 2 with margin, 
weight group 3 with margin, 
weight group 4 with margin, 
weight group 5 with margin, 
weight group 6 with margin, 
weight group 7 with msurgin. 



tons 
tons 
tons 
tons 
tons 
tons 
tons . 



I 

J 

( 
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3,15 Subroutine SEASPD 
3.15»1 Introduction 

This routine is based on the work of Michael Shen, 
Peter Gale and Kenneth Walker as reported in References 
4 and 5. In these references is given a computer program 
for computing the average sea speed. This computer 
program was altered by dropping the option to use a bow 
sonar dome and by deleting the read and write statements. 
The program was also rewritten in ASA standard FORTRAN. 
Except for these changes, subroutine SEASPD is the computer 
program developed by these three men. The interested 
reader is referred to these references for the development 
of the routine. 

Briefly, the routine assumes 1 (1) that the ship 
is operating in the North Atlantic Ocean over a long 
period of time 1 (2) that it is steaming into head seas 
one half of the time 1 (3) that the remainder of the time 
(when steaming at other headings) it is not forced to 
reduce speed due to ship motions. Part of the time the 
ship is steaming into head seas, its speed will bo limited 
be the sea state. In the routine this is simplified to ^ 
a single sea state limited speed. The remainder of the 
time the ship is assumed to be able to make its maximum 
sustained speed. The percentage of the time the ship 
is limited in speed is estimated based on average wave 
heights in the North Atlantic Ocean and on ship character- 
istics. 
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Tho average speed at sea is calculated assuming 
that the vessel will travel at its naximum sustained 
speed whenever it is not linited in speed by the sea 
state. When the vessel is steaaing at headings other 
than into head seas, its speed is not linited. This is 
assuned to bo one half tho tine. When steaming into 
head seas its speed will be linited only part of the 
time. Both the limiting speed and the porcentage of 
time that it applies are calculated in SEASPD. 

This routine has been included to indicate the 
effect that' changes in the vessels characteristics will 
have on its seakeeping ability. The many assumptions 
made in deriving the aver 2 ige sea speed make the speed 
predicted practically meaningless. However, the speed 
calculated can be used to compare two designs to decide 
which will give the better sea keeping performance. 

A flow chart for this routine is shown in Figure 46. 

3.15*2 Input List for Subroutine SEASPD 



B 




DPTRY 




CP 




H 




CWP 




LEH 




cx 




VS US 





5*15*3 Statement Descriptions 

For a detailed discussion of this routine see Refer 



encea 4 and 5* 



I 
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SEASPD SUBROUTINE 
FLOW -CHART 




I 
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Subroutine SEAS PD 

ij. . ^Ph 0 ;,pr^y. output is the value of average sea speed, 
AVSP. 

3.15*5 Nomenclature List 

The definitions used are the same as those given in 
References 4 eind 5 except for the followingi 



Value used . 
in SEASPD 


Value used 
references 


B 


BEAM 


LEN 


SL 


CWP " ' ' 


CW 


• 


B 


DPTRY' 


DISPO 


KY ■ ■ ' 


DKY 


lYY 


YYI 


f i'-i’rr.fA : " r ; tocv . . . 

JY 


YJ 


T t, .-.',>.vr»r l i 

lY 


YI 


vsus ■ 


WV 


vlr" 

> . ^ i 


VPX 


n ! A ‘v t, 




O]. 1 i M t 




;i*c- v, . .. . 




j p A* i J i * • - » i: « • . 




3 . i 6 . t If;; 





rou 
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3.16 Subroutine OUTPUT 

3.16.1 Introduction 

All program output is printed by subroutine OUTPUT 
except for error messages, A simplified flow chart for 
this routine is shown in Figtu^ 4?. Subroutine OUTPUT 
is divided into tv#o main parts. The first prints the 
input data and the second prints the output generated. 

If an error message is printed, only part of the 
output data is available and the data that is available 
is probably not correct. Therefore, when an error message 
occurs, only the input data is printed. Variable II has 
a value of 1 in this case. If 1;ho program successfully 
completes all proceeding steps, II is assigned a value 
of 2 and output data is also printed. 

Subroutine OUTPUT consists primarily of write and 
format statements. Some calculations are made but only 
to convert previously calculated values to a desired 
output form such as determining ratios between weight 
groups and full load weight. 

Only two decision points are used. One is the decision 
on whether or not to print output. The other checks the 
option used for inputting vehicle performance data and 
chooses the correct set of write and format statements 
to print the input data, 

3.16.2 Input for Subroutine OUTPUT 

All values which are printed are inputs to this 



routine , 
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OUTPUT SUBROUTINE 
FLOW CHART 




Figure 47 
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3«16.3 Statement Descriptions 

Detailed descriptions of the statements of this routine 
will not be given since most of the statements are self- 
explanatory. Use of the sample outputs given in Appendix 
B will be of use in understanding the format statements 
used. 

3.16.4 Output for Subroutine OUTPUT 

No values are generated in this routine for use 
elsewhere. 

3.16.5 Nomenclature 

All the variables used in this routine have been 
defined previously except for some dummy variables whose 
definition should bo clear from context. 



i 



i 



CHAPTBR IV 



EVALUATION OP THE PROGRAM 

A number of test runs were made on an IBM 370, model 
168, computer located in the Information Processing Center 
at M.I.T. The test runs vrere made using input data for 
four existing designs, the 378' WHEC, the 210' WHEC, the 
150’ WPB, and the Navy's Patrol Frigate design. The Cutter 
Model produces feasible designs for each of these ships 
except the 150' WPB. 

The design for the 150' WPB failed because a raised 
deck length greater than the ship length was required. 

The actual ship does not have a raised deck. The probable 
cause of this failure lies with the available area and 
volume relations used in the program. These relations 
were developed for ships with lengths greater than 400 
feet and probably have significant errors when extended 
to ships as short as 150 feet. Recommendations for cor- 
recting this error are given in Chapter V, 

A comparison between the results produced by the 
Cutter Model and the actual ships are shown in Tables 
6a amd 6b for the 378' WHBC and in Tables 7a and 7b for 
the 210' WMBC. The differences are reasonably small in 
most areas I however, a few do need discussion. 

The largest difference with both ships occurs in 
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Table 6a 

Comparison of Cutter Model Cutout to the Actual 
Ship for a 378* WHEC 







Beam 


Draft 


DO 


DIO 


D20 


Davg 


Disp. 




ft 


ft 


ft 


ft 


ft 


ft 


tons 


Ship 


350. 


42.0 


14.3 


34.75 


27.25 


27.5 


28,5^ 


2968.0 


Model 


350. 


43.32 


14.6 


38.33 


27.83 


28.18 


29.64 


3061.5 





WEIGHTS 


VCG 






Ship 


Model 


Ship 


Model 


WTGl 


917.46 


962.8 


19.12 


19.42 


WTG2 


383.6 


lU:t 


11.45 


11.52 


WTG3 


88.12 


23.49 


21.76 


WTG4 


80.2 


80.2 


31.56 


32.52 


WTG5 


268.8 


262.3 


19.43 


19.76 


WTG6 


325.1 


322.3 


25.63 


25.88 


WTG7 


39.0 


39.0 


35.26 


36.30 


WTLSHIP 


2102.2 


2215.0 


19.72 


19.98 


FUEL 


661.8 


653.2 


7.67 


8.35 


LUB OIL 


13.7 


16.0 


3.58 


3.62 


CREW 

PER EFFECTS 


;3 22.03 


13.7 

7.8 


^ 26.26 


27.00 

27.00 


PER STORES 


107.4 


95.1 


11.35 


14.36 


CARGO 

AMMO 


0. 

37.5 


0. 

37.5 


15.99 


16.42 


AIRCRAFT 


3.4 


3.4 


39.01 


40.08 


A/C FUEL 


19.6 


19.6 


9.56 


9.74 


WTFLLD 


2967.6 


3061.4 


16.61 


17.20 



Table 6b 
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Comparison of Cutter Model Output to the Actual 
Ship for a 378* WHEC 



Hull Volume 
Deckhouse Volume 

Total Arrangements Area 
Hull Arrangements Area 
Deckhouse Arrangements Area 

Office Space 
Messing Facilities 
Crew Special 
Officer S.R. 

Officer Sanitary 
CPO S.R. 

CPO Sanitary 

Crev/ Berthing 

Crew Sanitary 

CO S.R., Cabin & Pantry 

Commissary Stores 

Other Stores 

Workshops 

Repair Lockers 

Chain Lockers 

Uptakes 

Steering Gear & Windlass 
A/C & Fan Spaces 
I.C. Spaces 
Aux. Machinery Spaces 
Pilot house, Chartroom & CIC 
Passages 
Input Area 
•Excess Area 



Ship 


Model 




277600 


273734 


cu ft 


95611 


94705 


cu ft 


3O863 


29664 


sq ft 


19615 


I8522 


sq ft 


11248 


11142 


sq ft 


493 


510 


sq ft 


2713 


2713 


sq ft 


725 


787 


sq ft 


1157 


1168 


sq ft 


243 


192 


sq ft 


1040 


1040 


sq ft 


250 


208 


sq ft 


3641 


3168 


sq ft 


716 


720 


sq ft 


755 


753 


sq ft 


600 


652 


sq ft 


3015 


3231 


sq ft 


1161 


1254 


sq ft 


148 


100 


sq ft 


47 


40^ 


sq ft 


462 


0* 


sq ft 


736 


779 


sq ft 


369 


396 


sq ft 


312 


312 


sq ft 


i860 


2180 


sq ft 


1125 


1111 


sq ft 


2567 


2450 


sq ft 


5283€-» 


5900 


sq ft 


1445 


0 





♦ 460 sq ft included in input area 

*♦ Excess area was used for crew rec. rooms 
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Table 7a 
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Coiaparison of Cutter Modol Output to the Actual 









Ship 


for a 


210* WMEC 








ft 


Beam 


Draft 


DO 


DIO 


D20 


Davg 


Disp. 




ft 


ft 


ft 


ft 


ft 


ft 


tons 


Ship 


200. 


3^.00 


10.5 


20.75 


19.50 


20.90 


26.91 


990.0 


Model 


200. 


35.02 


8.98 


21.50 


19.50 


19.70 


24.97 


959.1 





WEIGHTS 


VCG 






Ship 


Model 


Ship 


Model 


WTCl 


379.5 


347.8 


16.4 


13.92 


WTG2 


102.6 


102.6 


9.5 


9.55 


WTG3 


^3.7 


75.4 


17.5 


22.28 


WTG4 


26.1 


26.1 


34.2 


33.31 


WTG5 


110.4 


109.6 


16.6 


19.88 


WTG6 


104.4 


106.8 


22.6 


25.93 


WTG7 


8.5 


8.5 


23.4 


23.21 


WTLSHIP 


775.2 


776.8 


17.1 


17.40 


FUEL 


152.3 


137.3 


8.17 


8.40 


LUB OIL 


3.8 


3.9 


2.5 


2.54 


CREW 


}8.6 


5.5 


20.43 


18.92 


FER EFFECTS 


3.2 


20.43 


18.92 


PER STORES 


47.3 


30.4 


10.99 


11.95 


CARGO 


0. 


0. 






AMMO 


2.6 


2.6 


8.4 


8.00 


AIRCRAFT 


0. 


0. 






A/C FUEL 


0. 


0. 






WTFLLD 


990.0 


959.5 


15.38 


15.87 






Ship 


Model 




Raised Deck Length 


164. 


120. 





Table 7b 



Comparison of Cutter Modol Cutout to the Actual 
Ship for a 210* W^C 



Hull Volume 
Deckhouse Vblume 

Total Arrangements Aree 
Hull Arrangements Area 
Deckhouse Arrangements Area 

Office Space 
Messing Facilities 
Crew Special 
Officer S.R. 

Officer Sanitary 
CPO S.R. 

CPO Sanitary 

Crew Berthing 

Crew Sanitary 

CO S.R., Cabin & Pantry 

Comaiesary Stores 

Other Stores 

Workohops 

Repair Lockers 

Chain Lockers 

Uptakes 

Steering Gear A Windlass 

A/C & Pan Spaces 

I.C. Spaces 

Aux. l!achinery Spaces 

Pilot house , Chartrooa & CIC 

Passages 

Input Area 

Excess Area 



Ship 


Model 






124700 


123049 


cu 


ft 


16100 


17840 


cu 


ft 


11242 


11021 


eq 


ft 


9345 


8922 


sq 


ft 


1897 


2099 


eq 


ft 


130 


192 


sq 


ft 


1452 


1381 


cq 


ft 


324 


313 


oq 


ft 


567 


657 


sq 


ft 


193 


108 


sq 


ft 


283 


320 


aq 


ft 


75 


64 


sq 


ft 


1658 


1254 


eq 


ft 


384 

406 


285 

398 


eq 

sq 


ft 

ft 


480 


393 


oq 


ft 


1294 


1233 


sq 


ft 


544 


382 


cq 


ft 


110 


100 


oq 


ft 


42 


40 


eq 


ft 


0 


0 






270 


390 


oq 


ft 


165 


151 


sq 


ft 


180 


180 


oq 


ft 


820 


832 


sq 


ft 


422 


486 


oq 


ft 


739 


910 


oq 


ft 


704 


704 


eq 


ft 


0 


247 


sq 


ft 
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weight group 3# electric plant. This difference results 
from the fact that greater generating capacity relative 
to the average load has been built into the model than is 
installed in existing ships. This increase in generating 
capacity reflects current Coast Guard policy for new 
designs and is not an error in the program. 

The differences in weight group 1 reflect the fact 
that the hull dimensions are not exactly the same for 
the program designs and the actual ships. Weight group 
1 is a function primarily of the hull cubic number. 

The fuel weight is calculated in the program based 
on specific fuel consumption, endurance range and endurance 
speed. An attempt has been made to approximate the actual 
fuel weight by varying one or more of these quantities 
but an exact match is difficult to obtain. 

The only other significant weight difference occurs 
in personnel stores which includes potable water. The 
model results are loss than the weights for the actual 
chips due to the fact that the ainoxmt of potable tTator 
per man has been reduced from the amount used on past 
designs. This is in accordance with current Coast Gutird 
policy for future designs. 

The vertical centers of gravity show quite a bit 
of variation when compared to the actual ships. However, 
the values in the program are only average values and 
will in general not agree with the values for any one 
ship. 
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In the case of the 210* WMBC, a raised deck is 
required. The model predicted a raised deck length 44 
feet shorter than that of the actual ship. This is another 
indication that the available strea and volume relations 
arc inaccurate for smaller ships. 

Hull areas are given in Tables 6b and 74. The areas 
used in the cutter model for berthing and sanitary spaces 
are based on current minimum Coast Guard standards. Theso 
areas in general are lo^r than the areas on current 
ships. 

Several runs were also made using data for the Kavy*s 
Patrol Frigate design. At the same time, runs were made on 
the Navy*s destroyer model, DD 07 , using similar data. 

The results of these rxins are shovm in Tables 8 and 9. 

The Patrol Frigate design is of interest because 
it is at the lower limit of applicability for the des- 
troyer model and at the upper limit of applicability for 
the cutter model. Agreement between the values calculated 
by the two models and the actual ship as shown in Table 
8, is reasonably good. Again, the group 3 weight is highest 
for the cutter model because 3 - 2000 KW generators are 
used instead of the 3 ” 1000 KW generators on the actual 
ship. Also, it appears that the greater emphasis placed 
on habitability by the Coast Guard is reflected in a 
greater group 6 weight. 

Table 9 gives the displacements predicted for greater 
lengths than the actual ship length. As the ship length 
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Comparison of Cutter Model Output to the Actual 
Ship and to DD07 Estimate for the PatroO. Frigate Design 





f^P 


Beam 


Draft 


DO 




ft 


ft 


ft 


Ship 


408. 


45.0 


14.5 




Model 


405. 


‘*5.95 


15.3^ 


41.07 


DDO7 


405. 


‘*6.72 


15.3 


41.00 



DIO 


D20 


Davg 


Disp. 


ft 


ft 


ft 


tons 








3517. 


28.92 


29.33 


31.01 


3706.1 


33.99 


34.80 


35.30 


3766.3 



WEIGHTS 





Ship 


WTGl 


1222.36 


WTG2 


256.1 


WTG3 


155.36 


WTG4 


92.58 


WTG5 


351.3 


WTG6 


268.82 


WTG7 


95.68 


WTLSHIP 


2442.20 


FUEL 


830.0 


LUB OIL 


20.22 


CPvEW & PE 


20.96 


PER STORES 


68.05 


CARGO 


67.0 


AJrlMO 


58.8 


AIRCRAFT 


10.0 


A/C FUEL 


0. 



3517. 



Model 


DD07 


1199.6 


1331.2 


258.0 


240.0 


211.2 


156.6 


70.7 


70.7 


317.5 


381.3 


333.8 


238.0 


85.7 


85.7 



2637.^ 2666.3 



818.7 


845.3 


20.2 


0.6 


20.9 


20.9 


66.6 


84.0 


75.0 


80.9 


58.8 


58.8 


10.0 


10.0 


0. 


0. 


3707.7 


3766.9 
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PATROL FRIGATE 
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Pull Load Displacement vs. Length Estimates 





Displacement 

Cutter 

DD07 Model 


405. 


3766.3 


3707.7 


415. 


3652.2 


3679.5 


425. 


3537.1 


3672.0 


435. 


. 3401.6 - 


3696.0 


445. 


3458.8 


3764.1 


455. 


3529.2 


3874.7 


465, 


3640.5 


3985.9 
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is increased, the displacements predicted by the cuttor 
model begin to diverge from those predicted by DD 07 . 

This indicates that the cuttor model begins to lose 
accuracy for ship lengths above 400 feet, the assumed 
upper limit of the cutter model. Table 9 also shows that 
the two models predict different optimums. This is a 
result of the assumptions built into the two models. 
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CHAPTER V 

RECOMMENDATIONS AND CONCLUSIONS 

For the 378* WHEC and 210* WMEC, the full load 
displacements predicted by the cutter model are within 
four percent of the actual full load displacements of 
these two vessels. This is bettor accuracy than is 
obtainable with other currently available models for the 
size range of Coast Guard vessels. However, only a limited 
number of designs were used to develop the cutter model 
and these same designs were used to check.it. The program 
should be used with caution until other designs are avail- 
able to confirm the relationships used in the model, 
t Although the cutter model appears to give accurate 
results in its current form, there are several areas 
where significant improvements could be made. The two 
most important improvements are i (1) develoiment of avail- 
able. arrangement area and volume relationships which are 
more applicable to ships of the size range used in the 
model. , and (2) expansion of the C^ array used in calculating 
horsepower, requirements. The electrical generator sizing 
criteria should also be given careful consideration to 
determine if »the sizes which result are what is really 
desired., - . . 

The-, develojaaent of area and volume relationships 



201 



J 

I 

i 

i 



202 

will require considerable developaental work. As a 
minimum, the variation of area with hiill depth, beam, 

Cp, Cx» tankage volume, and machinery box length must 
be considered. A computer program should be developed 
to calculate tankage volume and arrangements area for 
systematic changes in these variables. The results could 
then be converted into formulas for use in the cutter 
model. 

The expansion of the Cj. array could probably best 
be accomplished by fairing data from another model test 
series, such as series 64, into the data from Taylor's 
Standard Series, Where necessary some da^,.may have, to be 
assumed that is not available. An expansion of the 
array to speed length ratios of 3,0 and values up 
to 0.008 would be very desirable, Vfith the current 
array limits, many ships can not be run which may be 
very desirable. Also, if the sustained horsepower is 
given as an input, there are many times when the array 
limits are exceeded during an iteration even though the 
final speed value may be well within the limits of the 
array. 

Finally, the generator capacity predicted by the 
model appears to be excessive. This criteria should 
be given careful consideration. Existing 378' WHBC 
cutters seldom exceed the capacity of one 500 KW generator 
in normal steaming i however, the criteria built into the 
cutter model would result in 15OO KW generator sets being 
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installed. 

Some general recommendations on the use of the model 
are also in order. First, a ship synthesis model is only 
a design tool. Complete reliance should never he placed 
on a model I but, if chocked by hand calculations, the 
model will greatly inhance the engineer's ship design 
capability. Also, as a tool, the model must be kept 
current. This implies that periodic reviews must be made 
of the empirical relations used to insure that they are 
in agreement with current policy. Every new design built 
should also be plotted on the curves included in Chapter 
III as additional data. Changes should be made if the 
new data is not in agreement with the old. 

It is expected that the model will most often be 
used for designing new ships. Studied most often will 
probably bo changes in length and Cp. Those values have 
been arranged so that this can be done by changing only 
one input card. However, cheuigos in any of the input 
variables cein be investigated. 

If the model is used to investigate the sensitivity 
to changes in the design assumptions, it will be necessary 
to change statements in the program itself. If many of 
these studies are being performed, some of the fixed 
values in the program could be i^ead in on input cards. 

One constraint on all ship synthesis models is that 
variations in hull form can not be investigated. If 
hull forms much different from those of current designs are 
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desired* a nevf model has to he developed. It must be 
remembered that a ship synthesis model is really a multi- 
dimensional linearization of a complex problem. If 
large departures are made from the data base designs, 
the errors will increase in magnitude. 

Finally, the cutter model could be programmed as 
part of a ship optimization routine, but it is recommended 
that this not be done. Such a procedure puts the designer 
out of a controlling position while putting total reliance 
on an optimization criteria. Optimization criteria are 
difficult to define if more than just the cost of the 
vessel is in question. It is felt that ouch questions 
as the balance between payload and cost are best left 
to the designer and not to a rigid formula. 
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Description of Input 

Input for the cutter model can be divided into two 
main groups. The first group contains the data for the 
Cj. arrays used in estimating propulsion power. The second 
group contains the data for individual ships that are run. 

Data for the arrays is read only once regardless 
of how many test cases are run. This data precedes all 
other data and the values never change. The Cj. array 
data can be stored on tape or disk if desired. A list 
of the values used is given in Appendix D. 

The second group of data allows the user to specify 
the payload and performance as well as various other 
parameters of the ship the computer is to design. In 
order for the computer to produce reliable results, the 
input data must be the best that it is possible to specify. 
Past designs are a good starting point for determining 
reasonable input values. 

The minimum number of data cards for the second 
group of input is 14. These cards are described below. 
Where more than one card may be used this fact is noted. 

Card 1 — A single number is entered on this card 

in card column 1. The values that can be used 
are t 

0 - the program will stop 

1 - this is the normal value to use 



2 - a fine hull form will be assumed for 
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the first iteration, Cy = .0025. 

Card 2 — This card contains a heading that is 
printed above the output. All 80 columns 
can be used for any alphanumeric description 
desired. 

Card 3 — A single value is entered in card column 2. 
The total number of armament, aircraft and 
cargo payload data cards should be entered. 

A maximum of 20 can be used. 

Card 4 — This is really a set of cards containing 
the armament, aircraft £ind cargo payload data. 
The number of cards used is specified on card 3» 
Bach of these cards has the following format* 
(All data should contain a decimal point) 
columns 1 — 20 Any alphanumeric title 

describing the remainder of the data 
on the card. 

columns 21-25 Cargo weight, enter the 
weight in tons, 

columns 26-30 Cargo vcg, enter the ratio 
of cargo vcg divided by the depth of 
the ship at amidships, 
columns 31-35 Group 7 weight, .enter 



weight in tons 






«■: 
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columns 3^-40 Group 7 vcg, enter ratio 
of kg to depth at ajnidships. 
columns 41-45 Ammunition weight, enter 
weight in tons.. 

columns 46-50 Ammunition vcg, enter ratio 
of kg to depth at amidships, 
columns 51-55 Aircraft weight, enter 
weight in tons. 

columns 56-60 Aircraft vcg, enter ratio 
of kg to depth at amidships, 
columns 6l-65 Electrical load, enter value 
. in KW. 

columns 66-?0 Deckhouse Area, enter area 

required to he in deckhouse in sq, ft. 
columns 71-75 Hull Area, enter all remaining 
required area in sq, ft. 
columns 76 - 8 O Group 7 cost, enter acqui- 
sition cost in thousands of dollars. 

Except for the moment arms, the values given 
on these cards are summed by category and only this 
sum is used in calculations. Therefore it does not 
matter whether all the data is given on one card or 
on several. The only requirement is that the vcg 
values apply to the weights which immediately precede 
them. All values which are included in light ship 
weight group 7 must be given as data on these cards. 

If there is no payload data, card 3 should 
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contain a 1 and a blank card should be included 
for card 4, 

Card 5 — A single value is entered in card column 2. 
This is a count of the number of data cards 
described under card 6, A maximum of 20 can 
be used. 

Card 6 — This is really a set of electronics input 
data cards up to a maximum of 20, The number 
of cards is specified by card 5* 2ach of these 
cards has the following format t (All data should 
contain a decimal point) 

columns 1 — 20 Any alphanumeric title 
describing the remainder of the 
data on the card, 

columns 21-25 Group 4 weight, enter 
weight in tons, 

columns 26-30 Group 4 vcg, enter ratio 
of kg to depth at amidships, 
columns 31-35 Electrical load in KW. 
columns J6-U0 Deckhouse area, enter area 
required to be in deckhouse in sq, ft, 
columns 41-45 Hull area, enter remaining 
required area in sq. ft. 
columns 46-50 Group 4 cost, enter acqui- 
sition cost in thousands of dollars. 
Again the data may be broken up in any 
way desired as long as the correct sets of 
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weight and vcg axe given* •. On all of these 
caxds-and .on the cards .described under card 4, 
the-ar6a given ’.under 'deckhouse .area should 
included only .the . area .that , is required to be 
in the deckhouse* If an area can be either 
in the deckhouse or in the hull^ it should be 
listed under hull area. 

All items included in light ship weight 
group 4 must be entered on these cards. This 
includes standard navigation equipment. If 
there is no electronics data, card 5 should 
contain a 1 and a blank card should be included 
•w : for-card 6. 

Card 7 -- The number of officers, CPOs and crew 

members :are specified, on this card. All numbers 
must be' right- justified integers, 

'•a: •-'•jcolumns‘1-; — 5 'number of officers 
•'a: ' columns 6 — 10 number, of chief petty officers 
.-'^'•‘columns 11-15 number of enlisted men 
Card 8 — ' This • card’' contains cost" indices input. 

All values' should contain a decimal point, 

‘ • - columns' 1 5 labor rate in dollars/hour 

vr :-r •columns 6 — 10 group 1 index 
columns 11-15 - group 2 index 
' I : columns l6-20 group 3 index 
‘ « columns 21-25 group 5 index 
'••columns 26-30 group 6 index 



The values of the labor rate and cost indices 
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may be found under Section 3»1^» Figures and 

Card 9 — The number of endurance days for dry 

provisions and aircraft fuel payload are given 
on this card. All values should contain a 
decimal point. 

columns \ 5 endurance days for dry 

provisions 

columns 6 — 10 weight of aircraft fuel, tons 
columns 11-15 vcg of aircraft fuel, 

enter the ratio of kg to depth at 
amidships. 

Card 10 — A single value is entered in card column 

2. Two options are available. If a 1 is entered, 
use card 11a, This option allows the user 
to input sustained speed and have the computer 
calculate horsepower. If a 2 is entered use 
card 11 b. Here particulars of the machinery 
plant are input and the sustained speed is 
calculated. Any other value will have an 
unknown effect. 

Card 11a — Use this card type if a value of 1 is 
entered on card 10, The msiximum sustained 
speed, .endurance speed and range and propulsive 
coefficients at endurance speed and maximum 
speed are entered on this card. All data should 
contain a decimal point. 
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columns 1---5 
columns 6 — 10 
columns 11-20 
columns 21-25 
speed 

columns 26-30 
sppcd. 



sustained speed in knots 
endurance speed in knots 
endurance range , nautical miles 
propulsive coef. at endurance 

propulsive coef, at maximum 



Card 11b — Use this card type if a value of 2 is 

entered on card 10, The particulars of a specific 
machinery plant are entered on this card. All 
data should contain a decimal point. 

columns 1 5 maximum sustained horsepower, 

enter as HP/1000 

columns 6--10 minimum machinery box 
depth in feet 

columns 11-15 machinery box length in feet 
columns 16-20 total weight of group 2, tons 
columns 21-25 specific fuel consumption 
at max, horsepower in Ibs/HP-hr 
columns 26-30 specific fuel consumption 
at half power in Ibs/HP-hr 
columns 31-35 endurance speed in knots 
columns 36-45 endurance range, nautical miles 
columns 46-50 propulsion auxiliaries 
electrical load in KW 

columns 51-55 propulsive coef, at endurance 
speed 



columns 56-60 propulsive coef. at maximum 



speed 

columns 61-65 acquisition cost of group 2, 
thousands of dollars 
columns 66-70 lub oil weight in tons. 

Card 12 — This card contains various input options 
and dimensions. All values except the first 
should contain a decimal point. 

columns 1 5 MTYPE, a right- justified 

integer is entered here. 

If MTYPS = 0 the program returns 
to read a new card 1 
If MTYPE =1 a diesel plant is assumed 
if the machinery particulars have 
not been given 

If MTYPE =2 a CODOG plant is assumed 
If card 10 contains a value of 2, any 
value of MTYPS except 0 may be used, 
columns 6 — 10 free surface correction, 
enter value in feet. This doesnot 
have to be a positive number. A 
negative number can be interpreted 
as a KG below the average built into 
the program, 

columns 11-15 length between perpendiculars, 
enter the ships length in feet. 
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columns 16-20 Cp, enter a value for 
prismatic coefficient 
columns 21-25 enter a value for mid- 

ship section coefficient 
columns 26-30 GM/Beam, a value of 0.1 
is commonly used 

columns 31-35 design and builders margin, 
enter a value as a decimal, i.e. a 
10 % margin is entered as 0.1 
columns 35-45 this is the frictional 

resistance correction used in horse- 
power calculations. A value of 0.0004 
is commonly used. 

Card 12 may be repeated as many times as 
desired with new data each time. Each new card 
is interpreted as a new data case. 

Card 13 — If it is desired to input an entirely 

new set of payload and performance data or to 
stop, a value of 0 should be given in column 
5 of card 12. This will cause a new card 1 to 
be read. 

Card l4 — A new card 1 must be included. If new 
data will be given all of the data cards must 
be repeated, A value of 0 on card 1 will cause 
the program to stop. 
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A sample listing of the group 2 data cards is given 
below. The data used refers to a 210* W&ISC. 

Card 1 



Card 2 

210 WI.3C TEST RUN 



Card 3 
5 

Card 4 



Group 700 


4.7 1.64 






385. 76. 


Group 701 


.6 1,01 








Group 702 


2.8 .60 








Group 750 
Ammunition 


.4 .31 


2.56 

f 


.41 






f f 


t 


f t 


columns 


35 40 


45 


50 


75 80 



Card 5 

6 



Card 6 



Group 400 


.7 1 . 8521.99 319 . 190 . 


Group 401 


2.4 1.3 


Group 402 


.2 1.56 


Group 403 


10.4 .92 


Group 404 


9.1 2.69 


Group 450 


3.3 1.76 


lift t 

columns 25 30 35 40 50 


Card 7 


9 8 


57 


f 1 


t 


5 10 


15 ” columns 


Card 3 


4.75 1.4 


1.6 1.19 1.55 1.39 


f • 


• 1 1 1 


5 10 


15 20 25 30 - columns 



Card 9 
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60 . 0. 0. 

t f I 

5 10 15 “ columns 

Card 10 
2 

Card 11b 



4.5 19. 


32 .: 


102. 6 


.42 


.44 


14. 


4000. .71 


.55 .6 


599. 


3.85 


I t 


1 


t 


1 


t 


f 


t 1 1 


t t 


f 


t 


5 10 


15 


20 


25 


30 


35 


45 50 


55 60 


65 


70 


Card 12 




















1 - 1.61 


200. 


.59 


.892 


.1 


0. 


.0004 








t t 


f 


1 


1 


f 


t 


t 








5 10 


15 


20 


25 


30 


35 


45 - 


columns 







Card 13 

0 in column 5 
Card l4 



0 in column 1 
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CALL COST(LEN) XEC 0069 
CALL SEASPD(AVSP) XEC 0070 
CALL 0UTPUT(AVSPfEXCKG,2) XEC 0071 
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Group 1 " Hull Structure 



100 - Sholl Plating and Planking 

101 - Longitudinal and Transverse Framing 

102 - Inner Bottom 

103 - Platforms and Plats Below Lowermost Continuous Iteck 

104 - Fourth and Lower Continuous Decks 

105 - Third Dock 

106 - Second Deck 

107 - Main Deck and Hangar Deck 

108 - Forecastle and Poop Docks 

109 - Gallery Deck 

110 - Flight Deck, Landing Platforms, Special Purpose 

Itecks above Weather Dock 

111 - Superstructure 

112 - Foundations for Propulsion Plant Machinery 

113 - Foundations for Auxiliaries and Other Equipment 

114 - Structural Bulkheads 

115 - Trunks and Enclosures 

116 - Structural Sponsons 

117 - Armor 

118 - Aircraft Fuel Saddle Tank Structure 

119 - Structural Castings, Forgings and Equivalent Weldments 

120 - Sea Chests 

121 - Ballast and Buoyancy Units, Fixed or Fluid 

122 - Doors and Closures, Special Purpose 

123 - Doors, Hatches, Manholes and Scuttles - Nonballistic 

124 - Doors, Hatches, Manholes and Scuttles - Ballistic 

125 - Kingposts and Support Frames 

127 - Sonar Dome 

128 - Masts, Towers, Tetrapods, and Service Platforms 

150 - V/elding, Riveting, and Fastenings 

151 - Free Flooding Liquids 
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Group 2 - Propulsion 



200 - Boilers and Energy Converters (Non-Nuclear) 

201 - Propulsion Units 

202 - Main Condensers and Air Ejectors 

203 - Shafting, Bearings, and Propellers 

204 - Combustion Air Supply System 

205 “ Uptakes (Smoke Pipes) 

206 - Propulsion Contrcl Equipnent (Non-Nuclear) 

207 - Main Steam System 

208 - Feedwater and Condensate System 

209 - Circulating and Cooling Water Systems 

210 - Fuel Oil Service System 

211 - Lubricating Oil System 

212 - Nuclear Steam Generators 

213 - Reactors 

214 - Reactor Coolant System 

215 - Reactor Coolant Service Systems 
21 0 - Reactor Plant Auxiliary Systems 

217 - Nuclear Power Control and Instrumentation 

218 - Radiation Shielding (Primary) 

219 - Radiation Shielding (Secondary) 

250 - Propulsion Repair Parts 

251 - Propulsion Operating Fluids 



Group 3 - Electric Plant 



300 - Electric Power Generation 

301 - Power Distribution Sv/itchboairds 

302 - Power Distribution System (Cable) 

303 - Lighting System (Distribution and Fixtures) 

350 - Electric Plant Repair Parts 

351 - Electric Pov/er Generator Fluids 



Group 4 » Coiaaunication end Control 



400 - Navigation Equipmont ( Non-Elec tronic ) 

401 - Interior CoisEuni cation Systems and Equipment 

402 - Gun Fire Control System 

403 ” Countermeasure Systems (Hon -Electronic) 

404 - Electronic Countermeasure Syster^ (ECM) 

405 - Missile Fire Con+^ol Systems 

406 - ASM, Torpedo Fire Control Systems (Surface Ships) 
40? - Torpedo Fire Control System (Submarines) 

408 - Radar Systems 

409 - Radio Communication Systems 

410 - Electronic Navigation Systems 

411 - Space Vahiclo Electronic Tracking Systems 

412 - Sonar Systems 

413 - Electronic Tactical Data System 

450 - Communication and Control Repair Parts 

451 - Communication and Control Operating Fluids 
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Group 5 - Auxiliary Systeeis 



500 - Heating Systom 

501 - Ventilation System 

502 - Air Conditioning Systom 

503 - Refrigerating Spacos, Plantt and Equipment 

504 - Gasoline, HEAP, All Liquid Cargo Piping, Osygon® 

Nitrogen, and Aviation Lubricating Oil Systoma 

505 - Plumbing Installations 

506 - Firenain, Flushing, Sprinkler, Ifashdom, and Salt 

Water Service Systems 

507 » Fire Extinguishing Systems 

508 - Drainage, Ballast, I’rimmlng, Heeling, and Stabiliser 

Tank Syotons 

509 - Fresh V7ater System 

510 - Scuppers and Deck Drains 

511 - Fuel* and Diesol Oil Filling, Venting, Stowage and 

Transfer Systems 

512 - Tank Heating Systons 

513 - Compressed Air System 

51^ ~ Auxiliary Stoan, Exhaust Stoem, and Steam Drains 

515 - Buoyancy Control Systom 

516 - Miscellaneous Piping Systems 

517 - Distilling Plant 

518 - Steering Systems 

519 - Rudders 

520 - Mooring, Towing, Anchor and Aircraft Handling Systems 

and Deck Machinery 

521 - Elevators, Moving Stairways, Stores Strikodown, and 

Stores Handling Equipment 

522 - Operating Gear for Retracting and Elevating Units 

523 - Aircraft Elevators 

524 - Aircraft Arresting Goar, Barriers, and Barricades 

525 - Catapults and Jet Blast Deflectors 

527 - Diving Planes and Stabilizing Fins 

528 - Replenishment at Sea and Cargo H^dling 

550 - Auxiliary Systems Repair Parts 

551 - Auxiliary Systems Operating Fluids, Gases, and 

Stabilizer Fluids 
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Group 6 - Outfit and Furnishings 



600 - Hull Fittings 

601 - Boats, Boat Stowage, and Handling 

602 - Rigging and Canvas 

603 - Ladders and Gratings 

6 oh - Nonstructural Bulkheads and Nonstructural Doors 

605 - Painting 

606 - Deck Covering 

607 - Hull Insulation 

608 - Storerooms, Stowages and Lockers 

609 - Equipment for Utility Spaces 

610 - Equipment for Workshops, Laboratories, and Test Areas 

611 - Equipment for Galley, Pantry, Scullery and Commioeary 

612 - Furnishings for Living Spaces 

613 - Furnishings t Offices, Control Centers, Machinery Spaces 

614 - Furnishings for Medical, Dental, and Pharmaceutical 

615 - Radiation Shielding for Nuclear Support Facilities 

650 - Outfit and Fvimishings Repair Parts 

651 - Outfit and Furnishings Operating Fluids 



Group 7 » Armament 



700 - Guns and Gun Mounts 

701 - Anununition Handling Systems 

702 - Ammunition Stowage 

703 - Special Weapons, Handling and Stowage 

704 - Rocket and Missile Launching Devices (Surface to Air, 

Surface to Surface , and Sub-Surface ) 

705 - Rocket and Missile Launching Devices (Antisubmarine 

Warfare ) 

706 - Rocket, Missile, and Components Handling Systems 

707 - Rocket, Missile, and Components Stowage 

708 - Torpedo Tubes 

709 - Torpedo Handling and Stowage 

710 - Mine Handling Systems and Stowage 

711 - Small Arms and Pyrotechnic Stowage 

712 - Air-Launched V/eapons Handling Systems 

713 “ Air-Launched V/oapons Stowage 
720 - Cargo Munition Stowage 

745 - General Arrangement - Armament Drawings 

750 - Armament Repair Parts 

751 “ Armament Operating Fluids 
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